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An accurate QED calculation of transition probabilities for the low-lying two-electron configura- 
' tions of multicharged ions is presented. The calculation is performed for the nondegenerate states 

(ls2s)"^S'i, {ls2p3/2)^ Pi {Ml and A/2 transitions, respectively) and for the quasidegenerate states 
(ls2p) ^Pi, (ls2p) ^P\ {El transitions) decaying to the ground state (Isls) ^So- Two-electron ions 
^ I with nuclear charge numbers Z = 10 — 92 are considered. The line profile approach is employed for 

description of the process in multicharged ions within the framework of QED. 
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^ ■ I. INTRODUCTION 



Highly charged ions (HCI), in particular, two-electron HCI considered in the present work are under intensive 
experimental and theoretical investigation during the last decades. In HCI the electrons are propagating in the field 
of the nucleus, which exceeds in strength all other external electric fields accessible in laboratories. This allows for 
tests of QED in strong electric fields. The most precise experimental data have reached a relative level of accuracy of 



' about 1 % in one-electron ions (the measurement of the binding energy shift, i.e., the difference between the electron 

. ^ binding energy and the Dirac point-nucleus value for this energy, ground state, Z = 92 [T|]), about 0.4 % in two-electron 

^ ions (the measurement of binding energy of one electron in the two-electron ion, ground state, Z = 92 and about 

Oh. 



0.03 % in three-electron ions (the measurements of the 2pi/2 — 2si/2 energy difference for the first excited and ground 
states, Z = 92) Q. In the theoretical studies of the few-electron HCI such a level of accuracy requires the inclusion 
of the second-order (two-loop) radiative corrections as well as the screening of the first order (one-loop) radiative 
corrections. The nuclear size, nuclear recoil and even the nuclear polarization corrections appear to be of importance 
] as well. What concerns the interelectron interaction corrections, the first-, second- and partly third-order corrections 
QQ • should be accounted for in high-Z HCI. For intermediate Z values and, especially for the quasidegenerate energy 
\^ ' levels (see below), the inclusion of the interelectron interaction corrections to all-orders, at least within a simplified 
QO . treatment, becomes necessary. The existing experimental data for the transition probabilities are less accurate (3% 
■ for Z = 54 [^), but also require from theory to take into account the interelectron interaction and the lowest order 
, radiative corrections. Moreover, HCI can be used for the investigation of fundamental problems beyond QED: First, for 
' testing the Standard Model via the observation of Parity Nonconservation_(PNC) effects in HCI. Various suggestions 



OO ! on this subject were made in for one-electron HCI, in [1, 0, H, [^, [13, El for two-electron HCI and in for 
' four- and five-electron HCI, respectively. Regarding two-electron HCI such proposals are based exclusively on the 
crossings of the energy levels with opposite parity at some Z values. For example, according to the recent calculation 
[Hi the energy level splitting AE = E{2^Pq) - E{2^So) for He-like Gd {Z = 64), amounts to AE = 0.04 ± 0.74 eV. 
Since the PNC effect is proportional to = Re{A£: - iT/2}-^ = AE/{AE'^ + T^/i), where P = 0.016 eV is the 2^Pq 
^ , level width, the relative PNC effect can be unprecedently large (?/ = 0.05 for AE = r/2) or exactly zero for {AE — 0). 
To eliminate the uncertainty in calculations of the energy splitting AE the full account for the two-loop radiative 
corrections and the more accurate treatment of the interelectron interaction become indispensable. The evaluation of 
the PNC effects also demands a most precise knowledge of the transition probabilities (level width). Second, it was 
proposed to use HCI for the search of the variation of fundamental constan ts [Tsl. 1^ . Again the precise knowledge 



of the level crossings and transition probabilities is needed for this purpose [IS L 1X611. In many cases the two-electron 
ions are preferable for performing the corresponding experiments (see, e.g., Hal)- Therefore, it is necessary to 
develop adequate and accurate methods, which allow for predictions of energy levels, transition probabilities and 
other characteristics of two-electron HCI with utmost precision. 

The most investigated properties of HCI are the energy levels of the electron configurations. Since more than 40 
years in the numerous theoretical works a large variety of different methods based on the Relativistic Many-Body 
Theory (RMBT) and QED were suggested and employed in practical calculations. A short survey of these methods 
in a historical retrospective was presented recently in [17| . A common property of all these methods is that they are 
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exact to all-orders in the parameter aZ (a is the fine-structure constant and Z is the nuclear charge number). In 
many-electron systems the expansion in aZ implies an expansion with respect to the relativistic parameter vjc (v 
is the mean velocity of an atomic electron, c is the speed of light), so that the aZ-expansion methods can be only 
applied to nonrelativistic systems (low-Z atoms and ions). The Relativistic Dirac-Hartree-Fock (RDHF) method flSj 
and its natural extensions like Multi-Configurational RDHF (MC-RDHF) [l^ or a coupled-cluster method based on 
the RDHF approximation (CC-RDHF) [20|] were widely used in calculations performed for particular ion species. 
Within these methods the one-electron part of the many-body Hamiltonian is treated exactly as well as the Coulomb 
part of the interelectron interaction Hamiltonian. Only the Breit part of the interaction Hamiltonian is treated 
approximately. The validity of this approach was analyzed thoroughly in [U, [22|, [2^ . For few-electron ions with high- 
Z values the application of the perturbation theory with respect to the interelectron interaction becomes possible since 
the interelectron interaction is of the order of 1 jZ compared to the binding energy [23| . This feature of the many- 
electron atomic systems is exploited in the most powerful non-QED methods for the evaluation of the properties of 
HCI, in the Relativistic Many Body Perturbation Theory approach (RMBPT) By means of this method the most 
extensive calculations of the energy levels in two-electron HCI within a wide range of nuclear charge numbers Z were 
performed ^2^. Still, compared to the exact QED theory this method suffers from the lack of the negative energy 
contributions (though these contributions can be introduced with the special nontrivial efforts [27[ \A !29j). from 
the approximate treatment of the Breit interaction (without retardation) and from the neglect of the cross-photons 
interactions (which represent a special QED effect) in higher orders of perturbation theory. Moreover, the inclusion of 
radiative QED corrections within the RMBPT approach is possible only with the use of the aZ-expansion expressions 
(30| which, strictly speaking, are not valid for HCI. During the last few decades several rigorous QED approaches for 
the evaluation of the various properties of the HCI were formulated. Unlike the non-QED treatments, the application 
of QED allows for the consequent improvement of the accuracy of calculations. The first QED methods, based on the 
adiabatic S-matrix approach and the energy shift formula by Gell-Mann and Low [Slj (this formula was later adjusted 
by Sucher to the QED applications [32|]) was introduced in js^l and later applied to the various QED calculations 
by many authors. However, in higher orders of perturbation theory the adiabatic S-matrix approach becomes rather 
cumbersome due to the necessity of explicit evaluation of the adiabatic limit (when adiabatic parameter tends to 
zero). More advaiitageous for these purposes appeared to be the Two-Time Green's Function (TTGF) method first 
formulated in [s^ [Sa] (see also the recent review [3^). With this approach a large number of calculations concerning 
the higher-order (two-loop) radiative corrections to the energy levels 37] as well as the first-order radiative corrections 
to the hyperfine splittings in HCI [3^ and to the bound-electron g-factors in HCI 39] was per formed. An original 
approach with the covariant generalization of the evolution operator was recently developed in [40l |4l| . A special QED 
approach for the evaluation of the different characteristics of the HCI originates from the QED theory of the spectral 
line profile first developed by F. Low [13]. The application of the Line Profile Approach (LPA) to the evaluation of 
the energy levels shifts was first formulated in [43], simple examples were presented in [4^. The LPA possesses all 
the advantages of the other methods and allows for the evaluation of any higher-order corrections. A most general 
formulation of the LPA was given in [i^ |4^, [13] with application to the energy level calculations in HCI (see also the 
review [13 )■ In the present paper we apply the LPA to the high precision calculations of the transition probabilities 
in HCI. 

Because of inherent difhculties, the transition probabilities are less investigated than the energies. This can be 
explained by the presence of an extra photon line (emitted photon) in the corresponding Feynman graphs and by 
the poorer convergence of the QED perturbation theory. A calculation of the transition probabilities with respect to 
the relativistic corrections has been performed by Drake ^1, within the unified method. The work presents 
the first relativistic calculation of transition probabilities for the (ls2s) (Isls) ^5*0 transition. The work (ioj 

presents the calculation of transition probabilities for the (ls2p) ^Pi, ^Pi (Isls) transitions. A comprehensive 
review on the transition probabilities for two-electron ions has been presented by Johnson et al. 50] about one decade 
ago, where the transitions probabilities for low-lying two-electron configurations have been calculated for ions within 
the entire range of nuclear-charge numbers Z. In ISOl the RMBPT approach was employed. The contribution of the 
negative-energy states was discussed in '28', '29*1 . The first complete QED evaluation of transition probabilities in 
HCI with the account for the interelectron interaction and radiative corrections has been presented in (SI, , ,52] . The 
calculation was performed for nondegenerate states for transitions with emission of electric [5l| and magnetic [s^ 
photons, respectively. 

In this work we present a calculation of the transition probabilities for two-electron ions with nuclear charge 
Z — 10 — 100. The calculation is performed rigorously within the framework of QED. We also present a special 
technique developed to master the slow convergence of the QED perturbation theory in the case of the ions with 
intermediate Z values. The calculation is performed for the nondegenerate states (ls2s)^S'i, (ls2p3/2) '^-Pj {Ml and 
M2 transitions, respectively) and for the quasidegenerate levels (ls2p) ^Pi, (ls2p)'^Pi (£^1 transitions), decaying to 
the ground state (Isls) ^Sq. In the present work we apply the line profile approach (LPA) for the derivation of 
all necessary formulas and develop it for the evaluation of transition probabilities for quasidegenerate levels in the 
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framework of QED. In this paper we focus on the interelectron interaction corrections and leave the inclusion of the 
radiative corrections to subsequent studies. 

Our paper is organized as follows. In Section |ll] we present the general formulation of the LPA discussing its 
foundations and justification. A novel development of the LPA for the description of transition probabilities is 
presented in Section IIIII In Subsection IIII Al we consider transition probabilities for one-electron ions. The next 
subsections are devoted to two-electron ions. In Subsection llll BI the generic expressions for the transition probabilities 
in zeroth order (i.e., neglecting the interelectron interaction) are presented. The corresponding first-order expressions, 
where the one-photon exchange between the electrons is included, are given in Subsection IIII CI The formulas 
employed for the evaluation of amplitudes and transition probabilities are presented in Section [IVI The formulas for 
the nondegenerate case (Subsection IIV A[) and for the degenerate case (Subsection lIVBp are described separately. 
In Subsection IIV C] we derive the formulas for transition probabilities as they are applied in numerical calculations. 
Section |V] is devoted to the description of the computational methods. The discussion and analysis of the results, 
their comparison with the another available data and conclusions are found in the final Section IVII 



II. LINE PROFILE APPROACH 



The LPA is the version of QED perturbation theory (PT) which starts from the description of the atomic electrons 
as a set of noninteracting particles moving in the field of the nucleus y""'^ (Furry picture) and described by the 
solutions of the Dirac equation 

(TmP'^ - 70^^""' - m)7A = 0. (1) 

Here, are the components of the momentum 4-vector, p'^ — e is the bound electron energy, 7^ are the conventional 
Dirac matrices. Throughout this paper we use the relativistic units where H = c — 1 and the fine-structure constant 
a = /{he). The charge of the electron is e = — |e|. In this paper we designate the eigenvalues of Eq. ([1]) as £ while 
the physical one-electron energies as e = e + Ae. The idea of the LPA is to evaluate the corrections to energy (Ae) as 
the shift of position of resonance in some scattering process due to the interaction with the quantized electromagnetic 
field. This shift up to the very high orders of QED PT does not depend on the particular resonance process and when 
this dependence appears the concept of the energy level for the excited states cannot be strictly defined anymore 
[l3| • For the practical implementation of the LPA the process of the elastic photon scattering on atomic electron was 
employed. This procedure in the lowest QED PT order is depicted in Fig. [TJ 

According to the standard Feynman rules (see, e.i?., [23|). the S-matrix element for the graph depicted in Fig. [1] 
reads 

^('^ = {--ief j d''xud\di^ao{xu)l^-A^^j^>{x^)S{x^,xa)i^'^A^^^^^ (2) 

where = (t, r) denotes a spacetime point, ^aoi^) = V'oo ('")e~'^* is the one-electron wave function, ■i/J = ij;'^^'^ is the 
Dirac conjugated wave function and A^^"^\x) — A^^''^\r)e~^'^* is the 4-vector of the electromagnetic field potential 
(photon wave function), fc'' — (w, fe), A are the wave vector and polarization. The frequency of the absorbed and 
emitted photons are w = |fc| and lo' — respectively. We employ the standard covariant notations for 4- vectors 
together with the sign convention for the metric tensor {g^v) = diag(l, — 1, — 1, — 1). Einstein's sum convention is 
implied. The 4-dimensional volume is d^x = dtcfir. 

We employ the notations Xu, Xd for the "up" and "down" vertex coordinates in Fig.[T] These nonstandard notations 
will be convenient for the more complicated graphs considered below. The bound-electron propagator is represented 
in terms of an eigenmode decomposition with respect to one-electron eigenstates of Eq. ([T]) : 

2tt J ^ LJ„ - e„(l - lO) 

The sum over n runs over the entire Dirac spectrum. Note, that the subscript at the integration variable (a;„) is not 
the subject of summation over n. These indices are introduced again for the convenience in handling more complicated 
graphs. 

Inserting the expressions for the propagator and wave functions in Eq. ([2]) and integrating over time and frequency 
variables we arrive at 

5(2) ^ ^^27Ti)SiLu' -u;)e^y^ . (4) 
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Here we employed the shorthand notation 

= I d'rMr)rAi'-'\r)Mr)- (5) 



The amphtude {U) of the process of elastic photon scattering is related to the S-matrix element via 

S = {-27ri)S{Lu' ~ uj)U , (6) 

where uj and uj' represent the energies of the initial and final states, respectively. 

The resonance scattering means that the photon frequency lu is close to the energy difference between two atomic 
levels UJ — Ea — Sag + 0(0;), where a is the fine-structure constant. Accordingly, we have to retain only one term in 
the sum over n in Eq. Then the resonance amplitude looks like 



jj{2) _ ^2 ^aoa ^aap 



(7) 

UJ + Eao- 

This amplitude has a singularity &t lo — —Saa + £a- To avoid this singularity and to obtain the Lorcntz profile for 
the photon absorption and photon emission processes we have to consider the radiative insertions into the internal 
electron line in Fig. [1] This insertion is depicted in Fig. [51 
The corresponding matrix element can be written down as 



= {^lef j d''xud^xid^X2d^XdMxu)l^-A'^^l-'^'>{xu)S{xu.xi) 

xr'S{xuX2)l^^S{x2,Xd)r'A^^;^\xd)D^,^,{x^,X2)iJao{xd), (8) 
where D^^-^fj^^^xi, X2) denotes the photon propagator which in the Feynman gauge reads 

D^,^,{x,,X2) = ^/ rfr!/^,^,(|17|,ri2)e-'^"(*^-*^), (9) 

I^,^.A^,n2) - 6^''''^% (10) 

with ri2 = \ri — r2\ and the metric tensor g^^^2- 

Integration over the time and frequency variables and employment of Eq. ([6]) leads to the following expression for 
the amplitude in the resonance approximation 



UJ + Sao - £a 



(11) 



Here we introduced the energy-dependent matrix element of the electron self-energy 

Sw(0 = — / dn ^"""'^^'"'^ ^ (12) 

together with the shorthand notation 

IuiU2did2{^) = X! / ^''l^''2V'ui(T'l)'0«2(^2)7^'7'''-^MU'2(^>^12)V'di(ri)V'd2(''2)- (13) 



It is assumed that the ultraviolet divergent matrix element (I12p are renormalized in a standard way for the tightly 
bound electrons in atoms (see, e.g.^ [HI)- 

Repeating these insertions within the resonance approximation leads to the geometric progression. The resummation 
of this progression yields 

Ak',\')* Ak,X) 

U = e'^Jl^ (14) 

W-F£ao-K(w) 



5 



where 

KM - ^^w + AKH, (15) 

V^'^ = £a, (16) 
AVaiuj) = taaioJ + Sao) ■ (17) 

Within the resonance approximation we can replace uj + Sag by £a- Then AVa = Eaa(ea) represents the energy shift 
due to the lowest order electron self-energy correction. It is convenient to decompose explicitly the electron self-energy 
matrix element into real and imaginary parts 

Saa(ea) ^ La - '-T a , (18) 

where La is the electron self-energy correction to the energy of the excited level a and Ta is the width of this level. 
Thus, the pole in Eq. (|14p is shifted into the complex plane and the singularity on the real axis is avoided. The real 
parts of the resonance frequency in zeroth- and first-order of PT are given by 

^"■^^'("^ = -ea„+ea, (19) 

^rcs,(0+l) ^ ^Sao+ea + La. (20) 

Taking the amplitude Eq. ()14|) by square modulus, integrating over photon directions and summing over photon 
polarizations we arrive at the probability (cross-section) for the process of the photon scattering on an one-electron 
atom or ion. In the resonance approximation this probability factorizes into the product of absorption and emission 
probabilities with the same Lorentz-profilc factor [131 

'C('^'''=^) = 7 ^ft-TFJ' (21) 

where Taao is the partial width associated with the transition a — > oq. 

One of the most important problems within the LPA is the treatment of the ground state. An insertion of the 
electron self-energy correction into the outer electron lines in the standard Feynman graphs describing the elastic 
photon scattering on an atomic electron leads unavoidably to singularities. Since these outer lines correspond to 
the initial and final electron state, the problem of improving the energy of the ground state arises. Moreover the 
singularities leave the whole theory incomplete so far that the S-matrix for the bound electrons actually does not 
exist. One way to circumvent these difhculties was suggested decades ago by Sucher and Barbieri [53j whose idea was 
to evaluate the corrections to the transition probabilities via the imaginary part of the two-loop diagonal self-en ergy 
corrections for the ground state of an atom. Within the LPA we propose another solution of the problem [13, Hjj- 
To introduce the radiative corrections to the ground state we consider the more complicated two-photon process of 
the excitation of the resonant level a. This process starts from the one-photon absorption by an artificial "lower than 
ground" state a* (see Fig. [3]). This state plays the role of a regulator, which can be again removed at the end of 
evaluations. The insertions of the electron self-energy correction into the lower (or into the upper) electron propagator 
lead finally to the Lorentz line profile of the form [13, H^l : 



F +T 

(W - UJ^^^'Y + ^(Fa + F<j 



- aa) 



in the emission process a ^ ag -|- 7. In the most simple case of one channel decay the partial widths Faa^, Fa^a^ can 
be substituted by F^, Fa,,, respectively. In Eq. ((22l) the integration over the second photon frequency wq is performed 
and Fog represents the width of the level Cq. The presence of the width F^q is the only remnant of the introduced 
artificial "lower than ground" state a^,. In Eq. ([22|) uf^^ = —eag +ea includes both corrections La to the energy Sa and 
Lag to the energy e^o ■ If both states a and oq are excited states, the formula (|22p represents the known expression 
for the Lorentz line shape in case of the decaying final state. In case of the ground ag state, we may consider Tag as 
the regulating parameter. Setting Tag = at the end of the calculation, we obtain from Eq. ([22|) the emission line 
profile for the transition from the excited state a to the ground state ao with the radiative corrections for the ground 



state included in the definition of uj 



res 



Now being able to evaluate any desired property of an atom including the ground state energy corrections, we can 
disregard any insertions in the outer electron lines with their singularities and define the bound electron S-matrix in 
this way. Heuristically, this approach seems to be convincing, however, it looks unsatisfactory from the formal point 
of view. Therefore, we present a direct formal proof of the existence of the S-matrix for the bound electrons, based 
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on the adiabatic approach. This proof was given earher for the ^(oOjO) matrix [53|. However, ^(oOjO) matrix in 
QED provides difficulties with renormalization. Here we give the proof for S{oo, —oo) matrix for the first time. Our 
starting point is the Sucher adiabatic 5a„(oo,— cx)) matrix [sl] instead of Gell-Mann and Low adiabatic 5*^^(00,0) 
matrix [3l[, employed in [ssj. 

The standard description of an arbitrary process in the free-electron QED starts from the time evolution of an 
initial state to a final state governed by the evolution operator within the interaction representation based on the 
relation 

|$(oo)) = Texp|-i J dtH'^'^t)^ m^oo)) 

= S'(oo,-cx))|$(-oo)) . (23) 

Here, /f'"*(i) is the interaction Haniiltonian in the interaction representation, |"I>(±cx))) are the state vectors at 
asymptotic times t = ±oo and ^(oo, — oo) is the evolution operator, usually called S-matrix. The interaction between 
the particles involved is assumed to be absent a,tt = ±oo and the transition probabilities due to the particle interactions 
at the finite time moments can be expressed in terms of matrix elements of S-matrix. 

Contrary to this, the interaction between the bound electrons is permanently present. Within the adiabatic for- 
malism of Gell-Mann and Low [si'l the interaction Hamiltonian i?'"' (t) is replaced by the operator 

i7'"'(i) ^ e"^" 1*1 //'"*(<), (24) 

where Aq is the adiabatic parameter. Then, at the time moment t — ±oo the interaction is switched off and at i = is 
fully switched on. Using the interaction operator Eq. (|24p one can perform the QED calculations in the usual manner 
and then set Aa = at the end, thus restoring the full interaction for the entire time intervals. This allows for the 
extension of the established techniques for calculating free-electron S-matrix elements to bound electrons in atoms. 

Gell-Mann and Low [3l| derived a formula which yields the energy shift of bound electron states due to the 
interaction Eq. (I24|) in terms of the evolution operator S'a„(0,— oo). Sucher [33 | derived a symmetrized version of 
the energy shift formula, containing the matrix elements of the evolution operator S'a„(oo, — oo). On the basis of the 
Sucher's formula an adiabatic S-matrix approach for the evaluation of the energy corrections in bound electron QED 
was later developed [33j . 

Here we apply the adiabatic approach for another purpose. We will show that the singularities, arising after the 
electron self-energy insertions both in the initial and final outer electron lines in the S-matrix element corresponding 
to the Feynman graph Fig. [1] can be converted to the phase factor in the following way: 

Sx ^ {^2T:i)S{Lj' ~ Lu)Uexp I ) . (25) 

y iAa J 

This is an asymptotic equation (Aq —^ +0). The Aa-dependence is located in the imaginary exponent. Here the 
amplitude U differs from Eqs. (I6]),(l7]) by the replacement of Sao to — Sag + Sooao(£ao), where Saoao(£ao) is the 
diagonal matrix element of the lowest order electron self-energy operator. As for the ground state ag this matrix 
element is real, the phase factor does not contribute to the absolute value of amplitude defined by Eq. ((6]), and, 
accordingly, to the line profile. The proof is given in Appendix. This proof can be repeated for any QED correction of 
any order. This result justifies the employment of the energy (with the QED corrections included) instead of Sag 
in Eq. (jl4p . The remaining finite contributions from all the insertions in the outer lines present the QED corrections 
to the outer wave functions. Moreover, this proof can be repeated also for the few-electron atoms. In this case the 
ground state will be corrected not only for the QED corrections, but for the interelectron interaction as well. Thus, 
the problem of the existence of the bound electron S-matrix for the scattering process described by the Feynman 
graph Fig. [T] is solved in any order of QED perturbation theory. 

Now we go over to the most general formulation of the LPA. We will use the matrix formulation, that allows for the 
extension of the LPA to the case of quasidegenerate states. This formulation is also most suitable for the application 
of the LPA to the evaluation of transition probabilities. Within this formulation the generalization of Eq. (|14p looks 
like 

U = T+^T, (26) 

where the matrix T describes the absorption of the photon by the electron in the ground state oq with the excitation 
to the resonance (intermediate) state a. The matrix T'^ describes the emission of the photon with the transition 
a — > oq. The diagonal matrix (energy denominator) D{lu) is defined as 

Diu) = Lo + eag~~ , (27) 
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where m is the photon frequency. The resonance condition reads 

UJ'^' = -ta,+ea, (28) 

where e^o is the energy of the ground state, Sa is the Dirac energy of the state a. The energy is not necessarily 
equal to the Dirac energy, it may include already the radiative corrections (see discussion above). In Eq. (|27p the 
diagonal matrix D{u)) involves 

= Sa. (29) 

We employ here the matrix formulation of the single photon scattering amplitude (j26p on the one-electron ion in view 
of the further generalization of our approach to the quasidegenerate states in two-electron ions. This amplitude is 
described by the Feynman graph Fig. IH In this diagram we describe the photon interaction with the ground state by 
the boxes and deliberately omit the outer electron lines describing the ground state wave functions. The justification 
of this approach was given above. 

The next step is the insertion of the electron self energy corrections in the internal electron line within the resonance 



approximation (see Fig. [5]). The corresponding amplitude reads 47 1 



U = r+-— E(c. + £,„)-— T, (30) 
D{uj) ^[^) 

where S is the diagonal matrix corresponding to the regularized electron self-energy operator. In the case under 
consideration this matrix reduces to the diagonal matrix element of the electron self-energy operator for the state a. 

Continuing recursively this process, i.e., inserting two, three, etc, self-energies in the internal electron line in Fig. 
|4]and summing the geometric progression, we obtain finally 

where AV{uj) = Y.{ (jj ^ao)- Evaluation of the corresponding matrix element of the one-loop self-energy insertion at 
w = w"'' leads back to Eq. p8)) . 

Eq. ([3T|) illustrates the main idea of the LPA: the radiative corrections to the energy arise as the shifts of the 
resonances frequency due to the various insertions in the internal electron line in Fig.[3]in the resonance approximation. 
Graphically Eq. (|3ip can be represented by the Feynman graph Fig. [51 Instead of the correction E in the box in graph 
Fig. [6] any irreducible correction can be inserted; the corresponding energy shift will arise as the resonance frequency 
shift inj]q. J3ll). 

In [45i, |46|, the LPA was generalized to the few-electron ions, in particular for quasidegenerate states. The 
applications to two- and three-electron ions were presented. The general features of the LPA application to the 
evaluation of the transition probabilities were formulated in [l3| and exemplified with some numerical studies. 

In the next section of this paper we provide the detailed derivations of transition probabilities in the few-electron 
ions within the framework of the LPA. 



III. TRANSITION PROBABILITIES 



We are going to evaluate the transition probability for the process 

I^F, (32) 

where / is the initial two-electron state decaying to the final state F with emission of the photon luo. Within the 
framework of the LPA the state of an ion is associated with a position of the resonance. Therefore, we will consider 
a more general process which incorporates the transition ([5^ : 

A^^I^F^Ao, (33) 

i.e., a transition from the state Aq (let Aq be the ground state) to the state / with absorption of a photon to. Then, 
the state / decays to the state F with emission of the photon ujq and, finally, the state F decays back to the state Ao 
with emission of a photon lo' . The initial state (/) is associated with the resonance near u = —E^^ + Ef\ where Ef^ 
is the zero-order energy of the state / (sum of the Dirac energies). The final state (F) is defined by the resonance 
near w' = —Eaq + Ep\ The energy of the ground state Aq is given by Eao- 
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It will be shown below that in the resonance approximation the amplitude of the scattering process pS)) can be 
written as 

^ = Diu;') -AViu;') Die.) ^ AViu;) ^ ' ^^'^ 

The matrix T describes the absorption of the photon w by the ground state Aq, the matrix describes the emission 
of the photon uj' with the transition to the ground state Aq. The matrix D{lj) is defined by Eq. ([27|) . where F*^"^ 
is now the sum of the Dirac energies for the electrons which belong to the state /. The matrix D is diagonal in the 
basis of the two-electron functions in the j-j coupling scheme. The matrix of the interaction operator AV{uj) was 
investigated in [47| . Here we will construct it in the first order of the perturbation theory. 

The right denominator corresponds to the resonance associated with the state / and the left one defines the resonance 
for the state F. The function S(a;o) is a complicated vertex which describes the emission of photon wq by the ion in 
the state / decaying to the state F. The matrix element of the vertex ^(tJo) calculated on the eigenvectors <&/, 
of the matrices D{u}) — AV{uj) and Diuo') — AV{lo') corresponding to the states / and F, respectively, represents the 
amplitude of the decay process (|32| 

Ui^F = (S(wo))^^<j,^ . (35) 

The eigenvectors <l>i? and the vertex ^(wo) can be constructed order by order employing perturbation theory. This 
procedure is formulated consistently in following sections. 

Below in this section we will derive general formulas for the transition probabilities in two-electron ions in zeroth 
and first order of the QED perturbation theory considering the interelectron interaction as perturbation. 



A. One-electron ion 



In oder to introduce our notations we start from the one-photon transition in a one-electron ion within zeroth-order 
QED perturbation theory. This process with transition from the initial state / into the final state F is described by 
Eq. HMD. 

In zeroth-order QED perturbation theory the corresponding S-matrix element is given by Feynman graph depicted 
in Fig. [7] and reads 

S= [ d'^xi:F{r)e'*'''{-ie)Y'Ai^"'^"^*{r)e"^»'e-'*"^/ji{r). (36) 



The Dirac functions ?A/(t), '0f('") and the Dirac energies ei, Sf characterize the initial and final one-electron states. 
The emitted photon is described by the momentum 4- vector fco and the polarization Aq. 
In the coordinate representation the photon wave function 

^M(fc,A)(^)g-»^t ^ ./^e^(^)e-*('^*-'^'-) , = 1,2,3 (37) 

V LU 

describes a photon with the momentum k and polarization A (e*-"^^ is the polarization 4- vector). Here, the photon 
wave function is understood within the "transverse" gauge, which in the work [50.] is referred as the "velocity" gauge. 

Performing the integration over the time variable yields 

S = {-2TTi)6{eF+u;o-ei)e j d^r'4)F{r)-f^'A^i^°'^°>{r)il;i{r). (39) 

The expression for the amplitude (U) of the process is defined by Eq. Then, the amplitude corresponding to Eq. 
([39)1 reads 

U ^ e [ d^r^F{r)j''A^^"^^'>^*{r)i}i{r) . (40) 
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Within the framework of the line profile approach we consider a process described by Eq. (I33p . This process is 
depicted in Fig. [51 The S-matrix element corresponding to Fig. [5] is written as 



2,71 LUu^ ( 1 '^0) 

2tt ^ LJd - ed{l - iO) 

x(-ze)7''''4"/nrd)e-'"*''e-^*^(^»oVao(r-d) ■ (41) 

We employ notations "m", "c", "d" for the upper, central and lower vertices of Feynman graphs. Note, that subscripts 
at the integration variables w^, iOd refer to the corresponding vertices. After integration over the time variables and 
over the frequencies (a;„, ujd) we get the expression defining the amplitude {U) of the scattering process 



{-2T:i)5{uj' + ujo-u})U . (43) 



Being interested in the transition between / and F states we consider the frequencies of the absorbed and emitted 
photons satisfying the conditions 

w = -£ao+£/ + 0(a) (44) 

Oj' = -£a„ +£F + 0{a). (45) 

Let us also assume that the states /, F are well isolated. Hence, we can rewrite the amplitude in Eq. ([321) as 

U = / d^r.d^red^r,^,„(r„)7^"4t''''^*(^0 '^'^^l"^t^^''^^ 

X^M.^(^-o,Ao).(^^) V^/(^c)^/(^^) ^..^(/cA)(^^)^^^(^^)^^^ 
Eao + - £/ 

where R denotes the terms regular at w, uj' given by Eqs. ([331), ([3S1). The first term is singular and it defines the 
resonances corresponding to the initial state / and to the final state F . In the resonance approximation we retain 
only the terms singular at the positions of resonances what means that we neglect the terms denoted by R. The 
corresponding corrections are called nonresonant corrections. These corrections give rise to an asymmetry of the line 
profile and define the level of accuracy at which the concept of energy levels itself becomes inadequate for the analysis 
of experimental data. They are investigated in [s^ [s^l for highly charged ions and in [sl, [H, [60, [Ml [13, [H, [13] for 
the hydrogen atom. 

Aiming for the application of the LPA to two-electron ions we introduce the following notations. The vertex 
functions $ao, representing absorption of the photon by the electron in the state ao and emission of the photon 
with subsequent decay of an atom into the state oq, respectively, are 



Then, the expression for the amplitude takes the form 



<i>,„(r) = erAl^^^Hr)^Pao{r) (47) 
K,ir) = e^a„{r)rAi''-'>{r). (48) 



U 



e/d3,„,3,^d3r,<I>i(r-jM54£M 

J £aa + — £f 

X^.c^(/co,Ao)*(^^) V'/(rc)V^/(r-,) ^^^^^^^^^^ ^^g^ 
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Introducing notations 



Tnao ^ - j d>Vn(r)$ao(r), (50) 

Ti,n = - I d^r<^+{r)i^n{r) (51) 



we can reexpress Eq. (j46|) as 



1 



£ao +U)' - £F 



xe / d^r{pF{r)rAl^^«'^">{r)i^j{r) ■ Tja^+R- (52) 

J £ao + ^ — £l 

Introducing the matrices 

S„,„,(c^o) - e I d3rV;„,(r)7''4'=«'^°)*(r)V'„,(r), (53) 

D„,n2{^) = (^ +eao - £«i)^«i,«2 : (54) 

where (5ni,n2 means the Kronecker symbol, we can write the expression for the amphtude in matrix form 

U = T+p-\u;')E{u;o)D-\u;)Tao (55) 

" D[uj') ^{^) 

Expression (|56p coinside with Eq. ([M]) in zeroth order: AI/ = + O(q;). Taking into account the radiative corrections, 
the matrix AV^ will contain radiative insertions such as the self-energy and vacuum-polarization operators. As it was 
mentioned above in the present studies we will neglect the influence of the radiative corrections. 

According to Eq. ([55]) the amphtude of the process described by Eq. ((5^ reads in the zeroth order of the perturbation 
theory 

U = (S)^, . (57) 

As the matrix D{lu) is diagonal on the Dirac functions and AV^ = 0, the Dirac functions corresponding to the initial 
state / and to the final state F are the eigenvectors for the matrices D{uj) — AV{uj) and D{lo') — Ay(aj'), respectively. 

As a consequence of the application of the resonance approximation amplitude Eq. (|57p does not depend on the 
particular choice for the functions i i what means that the amplitude ([57)1 does not depend on how the initial 
state / was excited and how the final state F decayed. Accordinly, the state ao can be an arbitrary state. In the 
further derivations we will chose states $ao ; 'i'tio calculated within lowest orders of the perturbation theory. 



B. Two-electron ions: zeroth-order perturbation theory 

In the zeroth order of the QED perturbation theory the S-matrix element for the process Eq. (|32p in two-electron 
ions is given by the Feynman graphs Fig. [5) Within the framework of the line profile approach we consider process 
given by Eq. ([55]) . We assume the state as being the ground state. 

In the approximation of noninteracting electrons the S-matrix element corresponding to the scattering process (j33p 
is given by the Feynman graphs in Fig. [TOl The graphs (a) and (b) yield the same contribution, so we will consider 
twice the graph Fig. [TO] (a). 

In zeroth order the wave function of the ground state can be taken as Slater determinant 

<Hxi,X2) = i=det{V',„(a:i)V.o(x2)} = <^(ri,r2)e— (*^+*^). (58) 

Here, ipao{x) = i]jis+{x) = 'i/;is+(r)e~*'^^''*, iphQix) = ^pis-{x) refer to Dirac one-electron functions with different 
projections of the total one-electron angular momentum (identical with electron spin in case of Is-state). In zeroth 

order the ground state energy (Eaq) is = 2eis- The S-matrix element represented by the diagrams in Fig. [TOl is 
the same as the S-matrix element represented by Fig. [Sj and, accordingly, it is given by Eq. (|4ip . 
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With the purpose of employment of the matrix 5(wo) introduced in Eq. (|34p we rewrite the S-matrix element 
corresponding to the Feynman graph Fig. \TU\ (a) in the form 

S = ^^^^^ J d'^^nid^Xu2d^Xcid'^Xdid*Xd2duJuidLJdidLJn 
X _L \^ ''/'mi {'''ui )fpui (t'ci ) ^-iuj^, (t„, -t^, ) 



27r ^ - e„i(l - iO) 



ui 



27r ^ LUn- £„(1 - iO) 

xe~'''^^''-o+-^d{td, ~ idJ<^(r,,,r,J , (59) 

where we have introduced the two-electron vertex functions 

<^(n,r2) = e7^4'''H^i)*Ao(^i'^2), (60) 
<^(n,r2) - <)(ri,r2)e7^A;('=''^')(ri) (61) 

and function conjugated to the function ([58|) 



^^°„^(xi,X2) = ^det{V;ao(xi)^6o(x2)} = 'I'^°i(ri,r2)e--(*^+*^). (62) 

As in the case of the one-electron ion the functions <i>^^ , describe the properties of the scattering process (compare 
Eqs. (gaSHl)). 

The S-matrix element, corresponding to the Feynman graph Fig. [10] (a) contains two electron propagators. In Eq. 
([55|l we formally introduced the third propagator S{xu2,Xd2)- We have to show that the expression ([55]) leads to Eq. 
(|42| . With the aid of the Dirac-Sokhotsky formulas we can write 

[tj„ - e„(l - iO)]"^ — (5(cj„ - e„) - [-tj„ + e„ + zOe„)]"^ . (63) 

Employing this identity and integrating over variables 1^2 1 td2 ^-nd LOn, because of orthogonality of the Dirac functions, 
reduces the sum over n to terms n — Is only. The first term in the right-hand side of Eq. (j63p yields equation (j42p . 
while the second term vanishes after integration over ujui (because £„ > and the both poles lie in the same complex 
half plane). 

Going beyond the approximation of non-interacting electrons, the functions $Ao) ^Aq "wiH be more complicated 
functions than Slater determinants Eqs. (|60p . (j6ip . However, when calculating transition probability we can ignore 
the details involved in the preparation of the initial state (/) and the further decay of the final state (F) . Accordingly, 
in the resonance approximation, there is no need to specify the functions : ^Ao ■ 

In order to utilize efhciently the Feynman graphs technique within the framework of the line profile approach in the 
case of the two-electron ions we introduced two new elements: lower and upper boxes with letters Aq inside. These 
boxes describe the two-electron state Aq absorbing or emitting a photon. These lower and upper boxes correspond to 
the following expressions in the S-matrix elements (see Eq. ([5^ ) 

^-u,, iE.„ +-»5{td, - td2)'^Ao {Vd, ,rd2), 
{vu, , rn Je^*- (^-0 +-')5(t„^ - t„ J , 

respectively. Accordingly, in zeroth-order of the perturbation theory the S-matrix element for the scattering process 
Eq. ([55]) is represented by the graphs in Fig. [TT] and is given by Eq. ([55)1 . 
Integration over the time variables in (|59p yields the following expression 



S = i-z) 



(27r)^ ^ y dujui dujdi c^^»%5 



111 din 



'^'^Aouin^(^''''°'^°^*)iild.iTdinAo 

xK, -e„,(l-iO)]-i[L^d, ~erf,(l-iO)]-iK -£„(!- iO)]-i, (64) 
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which involves the following shorthand notations 

ZN[g5] = 5{Eao - i^ui - ujn)5{ujui - t^di + uJo)S{ujd^ + ujn - Eao - , (65) 
the complicated vertex 



Tn.n^Ao = -y ^-2 V'«i(n)V'n2(''2)$Ao(n,r'2) (66) 

and the one-electron matrix element 

Ai'^nl = Jd'ri.,,Ar)rA(,'''Hr)A^Ar). (67) 

We are interested in the transition probability at the frequencies corresponding to the positions of resonances 

^ = -Eao + Ef^ + 0(a) , (68) 
uj' = -EA,+Ef +0{a). (69) 

To zeroth order Ef'' = e/^ + e/^ and E^p'^ = epi + determine the positions of the resonances (sum of the 
Dirac energies) corresponding to the initial and final states, respectively. Within the framework of the resonance 
approximation one can retain only the terms which are singular at the positions of resonances. The Dirac energies 
E/j , S12, correspond to the positive-energy electron states, accordingly, within the resonance approximation 

we can omit all terms e„j < 0, e^i < 0, e„ < 0, what fixes the signs of the imaginary part of the poles in Eq. (j64p . 
Applying Eq. (|63p we can write the following equality 

^uju, - £„i (1 - «0)]-M^di - £di (1 - iO)]-i [^n - £n{l - iO)]-i 

= /N(^[£'ao - £ui -enT^[EAa + uj - Ed^ - £„] — -(5(£'ao + uj - ujd^ -e„) 

(70) 

The abbreviation Z ^ggj R for the product between Ag5] as given by Eq. ()65|) and the quantity R referring exclusively 
to the terms which are regular at the positions of the resonances is employed. These terms are regular, because the 
imaginary parts of the poles enter with equal signs. 

We employ Eq. ((70|) for the evaluation of Eq. ([64|) . Moreover, applying the resonance approximation in Eq. ((64|) 
all the regular terms R in (j70p can be omitted and the sum over ui,di,n reduces to the terms only, which satisfy the 
conditions 

Sd, + Sn = Ef^ = ei, + ei, , (71) 
Eu, + Sn = Ef = £Fi + ep^ ■ (72) 
Accordingly, for the contribution of the Feynman graph Fig. [TT] (a) we can write 
iS*' = (— 27ri)(5(aj — loq — lu') 

xT^lon, J^^o + ^' - en, - e„]-ie<°,f ")*[£;ao + ^ - - s^r^Td^nA, ■ (73) 

Here, we suppose that ui, di, n match with the conditions Eqs. (|7ip . (|72p . so that the index di runs over /i, I2, index 
Ml runs over Fi, F2 and index n runs over Ii, I2, Fi, F2. This yields a nonvanishing contribution only if /2 = F2 (single 
excitation) . 

In the same way the expression for S-matrix element corresponding to graph Fig. IllT b) can be derived 

S" = {—2'Ki)5{ijj — — i^') 

xrl„„„, [i?A„ + ^' - £n - + - e„ - ed^VT^^d^Ao , (74) 

where the states U2,d2,n now satisfy the conditions 

e„ + e<j, = Ef\ (75) 
+ = Ef. (76) 
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One can verify that the results for the graphs Fig. [TIT a') and Fig.[TTJb) are equal, i.e., 

= S\ (77) 



Our goal is to present expressions Eqs. (|73|). ([74|) in the form of Eq. ([34|) . In doing so, we consider the graph 
depicted on Fig. [121 The block S represents a complicated vertex describing the emission of the photon ujq. This 
vertex can be written in the form 

^{Xci, Xc2j , XS2) — ^{fcnfc2 7''^snfs2)^ S{tc2 ^ tci)S{ts-^ — tci)S{ts2 ^ tc^) ■ (78) 

The function S is a generic but yet unknown function. It can be derived under the requirement that the graph in Fig. 
[12] yields the same contribution as the graphs in Fig. [TT] The S-matrix corresponding to Fig. [12] appears as 

S = (-i)^ J d'^Xuid^Xu2d'^Xcid'^Xc2d'^Xsid'^Xs2d'^Xdid'^Xd2duJu^duJu2dL^didu}d2 
X $Ao (r-ni , r-„ Je**"i +"'^ <5(t„, - ) 

27r ^ - Em (1 - «0) 27r ^ - (1 - iO) 

^ ^ V'rfi(r,JV;rf,(rrfJ ^_,^^^ (,^^ -t.i ) A V ^'^2(^sjV^d.(rdJ ^„,^^^ (t^^ , 
27r 4^ ujd^ ~ Ed, (1 - iO) 2tt ^ ujd2 - £^2 (1 " iO) 

xe-**''i(^-*o+-)5(t,^ _t^J$^^(r.,^,rd,). (79) 

In the lowest order of perturbation theory the equations f75)) and (I74p follow from Eq. ([79)1 if we set 

S(r,,,r,,,r,,,r,J = 2e7^M|/;«'^«)*(r,JJ(r,, - r,JJ(r,, - r,J . (80) 

Consider now Eq. (I79p with S given by Eqs. ([75)1 and ([50)1 . respectively. After integration over the time variables 
we receive 

S = S^ + S' = Hf(^^^ ^^""^^ J ^'^"i^^^^'^^didwd^^ 

^-'AqMi«2 Midi "u2d2-^did2Ao 

x[uiui ~ £«i(l - iO)]"^[w„2 - £«2(1 - iO)]"^ 

x[cjd, -e<J,(l-^0)]-M^d2 -£d2(l-«0)]-S (81) 

where 

^^ (82] = '^(-E'ao + - t^ui - a;„2)'5(^"i + ^^«2 + ^^o - ^di - ^^2) 

X(5(tj<ji +t^<i2 - ^^Ao - t^) , (82) 

and ?iu2d2 is the Kronecker symbol. The employment of equalities analogous to Eq. (|63|) yields 

-£„i(l-iO)]"^K, -e„,(l-iO)]"^[wdi -£di(l -iO)]"^[wd, -£^,(1-^0)]"^ 

\^Aq + — £ui — £u2]^M-^Ao + ^ £di — £^2]^"^ (~^) '^(^«2 ~ £n2)'5(Wd2 — £^2) 

+A[g2]R. (83) 

The term is again understood as in Eq. ([70)1 above as shorthand notation for the regular part. 

Insertion oiEq. (|83)l into Eq. (|8T)) leads to the expression 

S" = (— 27ri)(5(a; — tJo ~ t^') 

^^Ao«iti2[-^'4o +t^'-£«i - £„2]~^2ui«2did2[-^Ao + W - - £d2]"^Tdid2Ao , (84) 

where 

^tiiti2dirf2 — 2eA|^^ 5u2d2 ■ (85) 



14 



We also suppose that E} ' = + e^a , E}, ' 

— £ui + £ji2 1 otherwise, this term is absent. Eq. (|84p together with Eq. 
© gives expression for the ampHtude of process Eq. ([55]) . With the use of Eq. ([55)1 one obtains the expression for 
the transition amphtude. In zeroth order the eigenfunctions are given by combinations of the Dirac functions 

in the j~j coupUng scheme. 

This was the goal of our derivations in this Section: to express the amplitude in the form equivalent to Eq. (j34[) . This 
presentation of the amplitude will help us to solve the problem of the transition probabilities for the quasidegenerate 
states. For the solution of this problem we need to present all the expressions in the generic matrix form Eq. ()34|) . 



C. Two-electron ion: first order perturbation theory (one-photon exchange) 

Now, we go over to the next order corrections to the transition probabilities and consider the one-photon exchange 
correction. This correction is represented by the graph in Fig. llSf a). The corresponding S-matrix element can be 
written as 



S = J d'^xid'^X2d'^Xuid'^Xu2d'^Xcid'^Xdid'^Xd2duJuj^cliOu2duJdidL^ 

„ J_ -ti) J_ ^U2{ru2)'>pu2ir2) i^^^ (t„^ _t,) 



27r ^ ujui - £«i (1 - iO) 27r ^ u)u2 - e«2 (1 " «0) 

1 27r cj^ £"7^(1 zO) 



27r 4^ cjdi - £di (1 - «0) 27r 4^ - (1 " «0) 

where ri2 = |ri — and the expressions for /^^^^ (|ri|, ri2) = ^^f^j ^'12) are defined in Coulomb gauge as 

= (87) 



or in Feynman gauge as 



ri2 axf ri2 

x(l - (5^io)(l - ^,.2o) 



/MiP2(f^,n2) = ^e^"'^-. (89) 
ri2 



We will also employ the following notation for the matrix element 

I%abi^) = y'd'nd3r2Via'(n)V;b'(r2)7r72^^/^fM2("'^i2)V'a(ri)V.(r2). (90) 

In Eq. (|86p again the additional electron propagator (sum over n) is artificially introduced with the same purpose as 
in the previous Subsection. 

Integration over the time variables in Eq. (|86p yields 

• \ 6 

I 



S = i-^f ( ^ ) (2^)' 



X J d^rid^r2d^ruid^ru2d^rcid^rd^d^rd2dL0uidLUu2dcudidLUd2duJndQ 
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X 



^ - (1 - iO) ^ - (1 - zO) 



Ml 

X 



^ - (1 - iO) ^ ujd2 - £d2 (1 - *0) 



where 



X(5(w„2 +f^-tJd2)^(-^^Ao - w + Wdi (92) 
Accordingly, in Eq. ([91]) we can set 

w„j = E'^n + - u;„2 , (93) 

ujdi = +L0 -L0d2, (94) 

uj„ = Eao +uj-uJu2- (95) 

Investigating the position of resonances near 

c^' = -EA„+Ef, (96) 

u = -Ea„ + Ef^ , (97) 

one can separate out those terms in Eq. (pij) that become singular near these resonances with the aid of the following 
sequence of equations (compare with Eq. (j63|) ) 



x[wdi - e<ii(l - *0)]"^[tJd2 - £^2(1 - «0)]"^ 

27r 

= I^^^Eao +w'-e«i - e„2]"^ — <^(^«2 -e"2)[-E^Ao + - -e«2]~^ (98) 
f 27r 

X { —S{uJd2 - £d2)[EAo +uj-ed,- ed2r^ 

-[Ea„ +uj~ujd2 - erfi(l - iO)]"^[-Wd2 +£d2 +iOed2r^ \ + 



Here R represents the terms which are regular in the vicinity of the resonances given by Eqs. (|96p and (I97|) . The first 
term in the curly brackets possesses a singularity (at the resonance Eq. ((97)) ) either in the case of Edi + Sd2 ~ Ef' or 
in the case of e„ + = Ej^'' . The second term becomes singular only in the latter case. 
Introducing the variable 

X ^ UJd2- £U2 (99) 

we rewrite the S-matrix element S^"^ corresponding to the Feynman graph Fig. [T3] (a) in the form 
S^'^ = (-27ri)5(w' + wo-w) 



X |e Inu2did2{\^d2 ~ |)[^Ao H" CJ — 6^ — ^ua] [^Aq -\- ^ — ^dx — ^^2] 
7> / ^7lU2 did2^'^\)\-^ ^d2 ~^ ^U2 ] 



2tt 



X [x- Eao -uj + edi+ £"2 - ^ [-E^a,, + w - £„ - £^2] ^ 



xTd,d2Ao ■ (100) 
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The first term in tiie curly brackets iias usually simple poles at two different points or it has a singularity of the second 
order if the points coincide. This term represents the first term of the geometric progression built for the initial state 
/ (see [l^l for details). Summation of the geometric progression results in a shift of the position of the resonance 
corresponding to the initial state and, accordingly, in a correction to the eigenvector of the initial state (^'/) (see Eq. 
psp ). As this term is taken into account while we generate the geometric progression, it does not contribute to the 
vertex operator and we can omit it here. 

Proceeding in a similar way for the evaluation of the S-matrix element S'''^, corresponding to the Feynman graph 
Fig. fTHT b). we get 



5"" 



(— 27r«)(5(aj' + loq ~ u) 



X [x-Eao -uj + ed2 - i^£d2V^[EAo + - -en]~^| 



(101) 



Note, that the graphs in Fig. fOfa) and Fig. \YMh) give equal contributions (for e^i + £^2 — ^ and 6^ + Sn = Ef \ 
respectively). Accordingly, the equality S^'^ — S"''^ holds. The vertexes corresponding to the graphs Fig. [TSKa) and 
(b) look like 



^(i)id ^ \^ pA^^^^^^^S 



e„+e^2.E\" 



-e^^ I dxInu2did2{M)[x - ed2 + £u2 - 'i-0ed2] ^ 



X [x- Eao - w + £di + £u2 - «Oedi] ^ 

)uiU2did2 7 



(102) 
(103) 



„(i)rd _ \^ eA'^''"-'^">5 

"uiU2did2 ^^^=^112" 



-e^^ J dx Iu^ndid2i\x\)[x ~ Edi + - iOedi] ^ 



X [x - Eao -uj + ed2 

C^(0)^(l)rdN 

V'—' JuiU2dia2 1 



(104) 
(105) 



where Sf^"^ is given by Eq. (f85|) . Here we have introduced matrices K'^^^^'^ and K'^^'^'^'^ for abbreviation of the matrix 
elements in the square brackets. Note, that if the initial state is well isolated the terms which do not match the 
condition e^i + £di = Ef^ are smaller by one order of perturbation theory and can be omitted. Accordingly, the 
matrices K become diagonal. 

As it was done for the zeroth-order corrections, the expressions (jlOOp and (jlOip should be cast into the form Eq. 
(p4)) . For this purpose consider the graph Fig. 14(a). The complicated vertex S has been already composed in the 
zeroth order (Eq. (|85p V Our goal is now to evaluate the interelectron interaction corrections to the vertex S. For that 
we will investigate the modification of a generic complicated vertex 'EP'^ after taking into account the interelectron 
interaction (Fig. 14(a)). The S-matrix element corresponding to Fig. 14(a) in the first order of perturbation theory 
in the interelectron interaction looks like 

S = (~*)^y d^xid*X2d^Xuid^Xu2d^Xcid'^Xc2d^Xsid^Xs2d^Xdid'^Xd2 

X dujui d^u2 duJsi duJs2 duJdi dwd2 ^Ag (^ui , f «2 ^^^^ ^"'-^ {tui - tu2 ) 



17 



27r ^ - (1 - lO) 27r ^ - (1 - iO) 

27r ^ t^,, - e,, (1 - iQ) 2ti ^ lUs, - Ss, (1 - iO) 

^ J_ V'di(ri)Vidi(r-dJ » V'rf2(^2)^d2(^d2) 

277 4^ - edi (1 - iO) 2tt ^ ujd^ - ea^ (1 - iQ) 

^e-~'i^^^^^o+^'^5{U, - U,)<^Ao{rd„rd,) . (106) 

Integration over the time variables yields 

■ \ 7 

I 



(''^'U2 )V'u2 (^^C2) ^gcn/„ „ „ „ N 

^ - £„i (1 - iO) ^ uju2 - (1 - «0) 

^ Wai - £si (1 ~ «0) ^ - £^2 (1 - «0) 

^ ^ f,Jj — — 7.11) -^^ — ^ /,Jj — I — 7, 1 



LJdl - £di (1 - «0) tJd2 - e<i2 (1 - iO) 



where 



S{fl + ujs2 ~ ujd2)S{^di + i^d2 - Eao - ^) • (108) 

Again we use the sequence of equations which separate out the terms {R) which are regular near the positions of the 
resonances under consideration (Eqs. ([96|l . ((97)) ): 

ATU5t^«i -e«i(l- jO)]"^K, -e„,(l-i0)]~^ 

x[cu,, - e,,{l - iO)]-^[LJs2 - es2{l ~ iO)]-^[u;d, - EdA'^ - iO)]-^[cJd2 - £d2{^ - ^OT^ 



= ^anill^T ) '^^'^"^ " £«2)'5(^S2 - £s2)S{oJd2 - £d2)5^^^^^^^ j,mS^^^^^^^^^mS^^^^^^^^^m 
x[Ea„ +uj' - e„, - £u2V^[Eao +i^-es,- es2]"M^Ao + w - e^i - £^2]"^ 

x[£;^„ + w - - es2]~^[EAo + ~ ujd2 - £rfi(l - »0)]"^(-l)[-tJd2 + £^2 +iO£d2]~^ 

+ ATMI ( ~ ) '^(^'"2 -£u2)5(wd2 -£d2)'5e„^+e„,,_Eg'>'^£rf,+£d,,_Ej">[^'4a + - -£„2]"^ 

x[-E^o +t^-t^s2 - £si(l - *0)]"^(-l)[-tJs2 +£s2 + jO£52]"^[£'Ao +^-£di -£^2]"^ 

Applying this result to Eq. (|107p and combining the graphs Figs. [T^l 14(a) (see Eq. (|84p ) for the S-matrix element 
S"^'*, corresponding to the Feynman graph Fig. 14(a) we get 

S"'^ = {-2TTi)S{uj' + cuo - uj) 
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~ l-l-^^gen 

'U2\ ^UlU2SlS2 



(l)siS2did2[^Ao + — ^di — £^2] ^ 

-e^^ J dx Is-^s^did2i\x\)[x - Sdi +£si -iO£di]"^ 

x[x~ Ea„ -uj + ed2+ £si " 'iO£d2]^^[EA„ +uj- Ss^- £52]"^ 



27r 



rfa;/sis2did2(kl)N - £di + £si + iOSsA ^ 



X[x~ Edi + Eao + - £s2 + «0£s2] ^ [Ea„ +UJ- Ed^- Ed2] ^ 



T, 



did^Ao 



The first term in the curly brackets has singularities at 

a;' = -Ea, 



-Eao + Edi + Ed2 



(110) 



(111) 

(112) 
(113) 



It can be considered as the first term of the geometric progression corresponding to the initial state. This progression 
can be summed up (see [47j). After this the position of the resonance corresponding to the initial state will include 
the interelectron interaction correction (one-photon exchange). 

The first term in the square brackets ((I)sis2ciid2 = <^sidi<^s2d2) represents the contribution of the graph Fig. [T^ 
The terms in the square brackets have singularities given by Eqs. (|111|) and by either (|112p or (|113[) . The last two 
terms represent the interelectron interaction correction to the generic vertex E!s°". The whole term in the square 
brackets corresponds both to the vertex S^''" and to the vertex Td^diAo^ which represents the process of excitation of 
the ground state by the photon w with transition to the excited state (/) . Accordingly, the contribution of this term 
to the vertex 2^°" appears with the power of 1/2. 

Suppose that the initial state is isolated, i.e., the admixture of the other states has a magnitude of the next order 
of the perturbation theory. Then, we can set 



Est + £5 



Edi + £d2 — E^ 



(0) 



(114) 



and omit the summation over si,S2 in the square brackets in Eq. (jllOp . Accordingly, the vertex ^'^'^ with the 
interelectron interaction correction given by the graph Fig. 14(a) will look like 



' 'U\U2d\d2 ' 'U\U2S\S2 



(l)siS2di(i2 



-e^— j dxls^s2dtd2{\x\)[^ ~ ^di + £si - iOEdj_]~^ 

X[x - Eao -UJ + Ed2+ £si - iO£d2]~^ 

-e^^ J dx Is-,s2did2{\x\)[x - Edi + ^si + «0£si]"^ 

X[x- Edi + Eao + W - £^2 + iO£s2]"^ 
l/2^ 

SiS2d-id2 



-,1/2 



(115) 



^U\U2SlS2 



cd 



The last equality defines the matrix K'^'^^'^'^. Note, that the correction factor appears under the square root. The term 
in square brackets in Eq. (jllOp can be equally referred to the vertex SS'^" and to the vertex Tdid2Ao- Accordingly, the 
part of this term connected with S^''", thus contributing to S'^'^, is the square root of this term. 
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-1/2 



The vertex E'^'^ should be equal to the sum of the contributions of Eqs. ([85)) . (|102l) . (|104p 

^cd ^ (^1 + + if(i)''*) . 

Accordingly, we derive (S'* S^''") 

= (^1 + X(i)''* + if (i^'^^j (^1 + /^(i^'^'^y 
Employing an expansion 

(l + x)-'/^ = l-^x + Oix^) (116) 
and neglecting the higher order terms we can write 

= (^1 + if (i)id + K^^^" - iif (i)-d^ + e 0(«2) . (117) 

Since we are interested only in the corrections of the zeroth and of the first order we can set 

uj - -Eao+Ei. (118) 

The first order correction to the S arise in the case of the reference state, i.e., when some of the following conditions 
are fulfilled 

Sn + £u2 = £tii + ^da ; (119) 

(120) 

£si +£s2 = £di +£d2 ■ (121) 

It is possible to describe the contributions of graphs Fig. [131 a). (b) as twice the contribution of the graph Fig. [TST a). 
Accordingly, the reference states are defined by Eqs. (|119p and (|12ip . respectively. 
Thus, we can write 

R:(l)ld , I^(l)rd _ 1 7^(l)cd ^ 2if (l^''^ - iif (^^'^'^ 

2 2 



-2e^^ J dxInu2did2{M)[x - £d2 + £u2 - iOi:d2] ^ 
'^^TT- / dx IsiS2did2{\x\)[x - +esi -iOedJ"^ 



dxlsis2did2{\x\)[x-edi+esi+i0esi] (122) 



This expression can be simplified with the aid of the following identities 

2Inu2did2{\x\)[x - ed2 + £u2 - «0]"^ 

-IstS2dxd2{\x\) {[2^-e<i2 +£u2 -«0]"^ + [a:-ed2 +£'«2 +«0]~^} 
= -fn«2di<i2(kl) {[2; " £^2 +e„2 -iO]"^ - [x~ed2 + £«2 +«0]"^} 

— j ^-?"««2did2(|2:|)(5(a; - £<j2 +£«2) • (123) 



Here the formula 



1 ^ '"%.) (124) 



(a; + «0)2 (a; — iO)^ z 9a; 

was utilized. 

The same procedure should be applied to Fig. [TST c'l.fd) and Fig. 14(b), where the one-photon exchange is inserted 
above the emission of the photon wq. 
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Finally, we can write the following expression for the vertex S 

5 = S(«)+S(i)+eO(a'), (125) 



where 



7(0) 

~'uiU2did2 

7(1) 

~'uiU2dld2 



V g3 4(feo,Ao)* CI IN 



6 ^^-^MlM2nd2 (I'^l) 



(fco.Ao)* 
ndi 



(126) 



(127) 



2-»i»iTt.U2 



Eq. (|127p represents the reducible part of the first-order corrections, i.e., the reference states contribution. 

Having constructed a general expression for the vertex S we can apply the formulas derived for the generic graphs 
Figs. 14 and[T2lfor evaluating the contributions of the graphs Figs. [Til and fT3l Now we can express these contributions 
via the matrix S, what enables us to extend the calculations to the quasidegenerate levels. To derive the formula 
for the amplitude Eq. p5)) . we will have to consider separately the case of nondegenerate levels and the case of 
quasidegenerate levels, respectively. 



IV. EVALUATION OF TRANSITION PROBABILITIES 



Evaluating the transition probabilities we should distinguish nondegenerate and quasidegenerate levels. For the 
nondegenerate levels standard QED perturbation theory can be applied. Configurations are called quasidegenerate, 
if they can not be considered as being well isolated. For these configurations the interelectron interaction must be 
taken into account up to higher orders. Accordingly, this requires to develop a special technique. 

In the next two subsections we derive expression for the amplitude of the scattering process for nondegenerate 
and quasidegenerate levels, respectively. In the last subsection we write down the final expressions for transition 
probabilities suitable for numerical calculations. 



A. Nondegenerate levels 



Here, we will suppose that the initial and the final states are well isolated. The set of graphs in Figs. [TT] and [T3] 
should be divided into two subsets: reducible (containing the reference states) and irreducible. For the zero-order and 
the reducible subset of diagrams formulas Eqs. ([55]) . (|125p can be applied, where the functions $/, <I>i? are given by 
a combination of the two electron determinants in j-j coupling scheme. For irreducible subsets of diagrams we can 
apply a procedure described below. 

Consider the first terms in the curly brackets in Eqs. pOOp and (jlOip 

5'^ = {-2m)S{cj' +u;o- c^)r+ [Eao + ^' - £u, - e^.^' (128) 
^ + ^ - - ^n]-'lnu2d^d2{\ed2 - £«2l)| 

X [Eao + ^ — £di — £d2]^^Td^d2Ao J 



-U2 \ 



x[EAo+u)-ed,-ed2r^Td,d2Ao- (129) 

The prime at the sum symbol indicates that in Eq. (jlOOp terms for which Sd^ + Sd2 — £n ~ £u2 = (and in Eq. (|10ip 
terms, where Sd-^ + £d2 ~ £«! ^ = 0) holds are omitted. Since the levels /, F are well isolated the expressions in the 
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curly brackets in (|128p . (|129p can be viewed as corrections to the vertex S; one can also set lu = —Eao + Ef''^ in the 
vertex. Thus, we have to take into account the following corrections 

^ — ( 1)6 ^ ' ^liri \^di ~^ -fnU2(ii!i2 (I^<i2 ^«2l) 

^di +'d2 -»ii-eu2^'' 

+ (-1)6^ Al^,"/"^*[edi +ed2 -en]"^^«in<iid2(k«i -^dj) (130) 

(this is the contribution of the graphs Fig. [T3Ka),(b)) and 

^(l)U _ \^ J rie- _ e- lUc- -L e _ e- _ c- 1 -1 /I C^o , Ao)* 

" — ^l'-- / J ^Iiiii2nd2\\^d2 ^U2\)i^ui "T tdj tnj -^din 

eiil+eu2"^""=<i2'^° 

+(-l)e3 /u,„2d,„(ku, -edJ)[eu,+£u2-ed,-e«]~Mi'°/"^* (131) 

(this is the contribution of the graphs in Fig.[T3Kc),(d)). Accordingly, for nondegenerate levels the amplitude Eq. ([55)1 
is given by the matrix element of 

S = s(")+S«+sWd + 5(i)u (132) 

evaluated with the aid of the zeroth-order wave functions corresponding to the states / and i.e., by means of 
two-electron determinants in the j-j coupling scheme. 

B. Quasidegenerate levels 

In the previous section we introduced the vertex S via the expression Eq. (|34p . In order to derive the amplitude 
as defined in Eq. ([55]) the wave functions <&/, have to be constructed. These functions are eigenvectors of the 
matrix V which was investigated in [47| . Diagonalization of the matrix is a serious task, because V has infinite 
dimension. One possible solution of this problem is the substitution by a large but finite matrix. Another strategy 
is the modification of a perturbation theory. Here we will concentrate on the development of a proper perturbation 
theory. 

The perturbation theory for the case of a nondegenerate level (as well as for the case of the fully degenerate levels) 
is well known [6^1 • Here we will apply it to the case of the quasidegenerate levels. Considering N two-electron states 
{^} defined in the j-j coupling scheme we assume that these states are mixing with each other, i.e., they have 
the same symmetry and the corresponding energy levels are close to each other. Under such condition the standard 
perturbation theory may not work and we have to modify it. These N states compose a set g = {'^i^,ig = 1, . . . , N}. 
The idea is now to build an eigenvector corresponding to a state ^„ £ g. We also suppose, that all the other 
states (beyond the set g) are either nonmixing with the state Ug or their energy levels are far enough from the level 
Ug, i.e., that the set g is large enough to incorporate all the closely lying levels. Then perturbation theory will again 
work. Otherwise the set g has to be enlarged. Similar, but not equivalent schemes were considered earlier in the 
frames of RMBPT (6^. Here, we apply it for the first time in QED. 

It is convenient to write the matrix V in a block form 



V = 



V21 V22 



(133) 



where the block Vn is constructed entirely on the states from the set g and the block V22 does not contain states 
from the set g. The matrix V can be decomposed as 

V = V'-^^ +AV, (134) 

where V^^^ is a diagonal matrix (a sum of the Dirac energies). The matrix AV contains the small parameter a (the 
expansion parameter of the QED perturbation theory) and can be treated as a perturbation. In what follows we 
restrict ourselves to the interelectron interaction corrections. In the lowest order these corrections reduce to the one 
photon exchange correction 

AV = Y I^i\b-b'\)a'b'ab, (135) 
g=c,t 



22 



We can write the matrix V as 



V = 



>11 


V12' 




V21 


V22 





y(0) 



^V2l 



V, 



(0) 



12 

AV2 



22 



(136) 



The block matrix Vn is finite and can be diagonahzed numerically according to 

B+ViiB. 



T/diag 



(137) 

Since in general is a complex valued symmetrical matrix, i.e., Vij = Vji the matrix i? is a complex orthogonal 
matrix 



B*B = /. 

Here / is a unit matrix [lij = 5ij) of the proper dimension. The superscript t in Eq. p38|) means transposition. 
Compose a matrix 



which is also an orthogonal matrix 
Acting by the matrix A onV yields 



A = 



A' A 



B 
/ 



V = AWA^ 



V^^^^ B'AVi2 
AV21B V2 
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(138) 

I. 

(139) 
(140) 

(141) 



Since we have supposed that the required state Ug is weakly mixing with the states not included in the set g, the 
matrix V can be diagonahzed with the standard procedure [65j 

ydiag ^ C'VC, (142) 

where the matrix C can be built order by order. The zeroth and the first orders of the matrix C look like 

(B'Ayi2)., 





(Ay2i-B)i 



Ej-Ei 
E,-Ei 



The diagonahzed matrices V and V coincide, so we can write 

ydiag = |/diag ^ ^ACfV{AC) . 
Accordingly, an eigenvector $ corresponding to a basic function ^' can be defined as 

$ = AC^ . 

Now we represent the state rig e 5 in terms of a perturbation expansion 



AC-^r. 



kg eg 



(0) 



B, 



21 )kl„ 



(0) 



E 



(0) 



(0) 



(143) 

(144) 
(145) 

(146) 



An expression for AV21 is given by Eq. p35p . Summation over index k means the summation over all two-electron 
configurations coupling scheme) including the negative part of the Dirac spectrum (not included in the set g). 
The employment of the j~j coupling scheme is not obligatory here. 

In case when the investigated state Ug is well isolated nondegenerate level and the set g consists only of this single 
state, i.e., g — {^'ng} the matrix B is just a one-dimensional unit matrix. It is easy to ensure that formula (I146|) 
together with Eq. (|125p gives the same result as Eq. (|132p (taking into account only the zeroth- and first-order 
corrections) . 

We again would like to point out the QED effects that are now taken into account in the framework of the LPA, 
and which are missing in the Rclativistic Many Body Perturbation Theory (RMBPT) [s^]: The first is the inclusion 
of the retardation (see Eq. (|135p ). the second is the account of the negative part of the Dirac spectrum (summation 
over k in Eq. (|146p and over n in Eqs. (|130p . (|13ip ) and the third is the incorporation of a nonzero contribution of 
S^-'^) in Eq. p26p (reference state contribution). 

The amplitude (U) of the scattering process for quasidegenerate levels is given by Eq. ([55)1 where the eigenvectors 
are defined by Eq. (I146P and the vertex operator is given by Eqs. (|126p . (|127p . 
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C. Transition probability 

Based on the scattering amplitude U of the process I ^ F with emission of the photon lo the transition probability 
between / and F states is given by formula 



(147) 



where u = k/\k\. Ej, Ep are the energies of the initial and final states, respectively. These energies comprise the Dirac 
energies and the one-photon exchange corrections. For quasidegenerate levels they are given by the corresponding 
eigenvalues of the matrix V (Eq. (|134p '). The photon frequency a; = |fc| should be set equal to w = _E/ — Ep. Eq. 
(I147p defines the full transition probability, i.e., integration over all momenta of the photon (fc) and summation over 
all polarizations of the photon (A) is performed. 

The integration over k and summation over A are taken analytically. In Eq. (|147p only the photon wave functions 
depend on k and A. Accordingly, in a very general way we can consider the case when 

C/(°) = A^^;i;> (148) 
d'rV;„,(r)7''4''^^*(^)V^".(^)' (149) 



where A^^'^\r) is given by Eqs. ([38]) . (l37|) . The corresponding expression for the transition probability can be written 
as [63 



(2^)2 



{AZ^*{r,u)U^,\ . (150) 



Here, the notation introduced in Eq. (|67p is employed. The summations run over the angular momenta of the photon 

(j) and projections (to). The 4-vector ^(*^'^)'^ = (V, A) corresponds to magnetic (Af) and electric {E) photons, 
respectively. In the case of the magnetic photons 

V^^\r,u:) = 0, (151) 

Af^\r,u) = ^g.^^^r)Y,,„,{n). (152) 

With the appropriate choice of the gauge for the electric photons we can write 

Vl^\r,u) = (153) 



In Eqs. (fT^ . pM)) the radial functions 

9i{x) = 4^\/^^'+i/2(2^) (155) 

involve Bessel functions Ji+i/2{x) is the of the first kind [68], Yjim {I = j — + 1) denotes the vector spherical 

harmonics [13, '69j depending on angles n = r/\r\. Formulas (jl53l) . (|154|) correspond to the photon wave function 
given by Eqs. ([38]), ([37]) . specified within the "transverse" gauge ^50] . 

For the nonrelativistic limit the more convenient gauge is represented by the transformation A ^ A + vxik^t), 
V ^ V + x{k,t) with 



X(fe,t) = (5(c.-|fe|)^^y,„,(i.)e-^"*, (156) 

where Yj^ii') is the spherical harmonics ^6911. This transformation affects only electric photons. Accordingly, in the 
nontransverse gauge the 4-vector A'^-^^ appears as 



vl2\r,Lu) ^ ^^^Jl±lg^ic,r)Y,M, (157) 
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27r /2j + l 



5j+i(wr)y,j+i,„(n) , 



(158) 



In the work [50] this gauge is referred as "length" gauge. 

Comparing Eq. (I147P and Eq. p50p we can express the transition probabiUty in terms of the corresponding scattering 



amphtudes t^j™'*^^ as 



jrn ^ 



(159) 



where A^'^''^"' are substituted by A^^l'^"^ , respectively. This expression was applied for the numerical calculations of 
the transition probabilities. 

The modified amplitudes J/j^'*^'' are derived within perturbation theory. As we take into account only the correc- 
tions up to zeroth and first order, i.e., 



U. 



(E,M) 



then, the squared absolute values of U. 

2 



(EM) 
jm 



jAE.Mm jj{E.M)(l) 



u. 



{E,M) 



U, 



[EMm 

jm 



read 

2 



"} 



{E,M){1) 



(160) 



(161) 



The last term in Eq. (jl6ip already corresponds to a corrections of second order and can be, accordingly, disregarded 
in the calculations. However, this term may serve as an estimate for the magnitude of the higher-orders corrections 
(i.e., for the error magnitude). The contributions of this term are given in the tables as AM^(2+' . 



V. NUMERICAL METHODS 



In the numerical calculations an ion is considered to be enclosed into a spherical box with the radius R = 60/ (aZ) 
(in the relativistic units), where a is the fine-structure constant, Z is the nuclear charge. The size of the box reflects 
the size of the volume, where the physical processes of the interest (photon emission, interelectron interaction) mainly 
occur for the two-electron ions with high and intermediate Z. Hence, the electron spectrum becomes discrete. The 
Dirac spectrum in the external field of the nucleus is constructed in terms of B-splines [tO, [Zli ■ We used B-splines of 
order 8 and a grid with 50 knots. 

Expression Eq. p46p for the eigenvectors <!>„ involves the zeroth- and first-order terms of the perturbation expan- 
sion. The matrix B required for the calculation of the eigenvectors $„g was generated perturbatively. For a given set 
(7, these two perturbation series are independent. 

The matrix V employed in the construction of the matrix B was borrowed from our work 47] , where it was evaluated 
up to the second order with respect to the interelectron interaction corrections. The matrix AV21 involved in Eq. (|146p 
is given by Eq. (|135p and includes only the first order of the interelectron interaction corrections. It was calculated 
in the present work. 

The spatial integration in the matrix elements of the type Eqs. ([S7|) and (|90p is performed in spherical coordinates. 
The integration over the angular variables can be performed analytically, while the integration over the radial variables 
is performed numerically. For the numerical integration we employed Gauss-Legendre quadratures, which yield a 
numerical accuracy of our calculations about 0.03%. 



VI. NUMERICAL RESULTS AND DISCUSSION 



In Tables Ill lIVI we present numerical results for Ml, M2 and £^l-transition probabilities for low-lying two-electron 
configurations in HCI. The values are given in units s^^ and the digits in square brackets refer to the power of 10. 

In Tables HI and HIl we consider transition probabilities between the (ls2s) '^Si configuration and the ground (Isls) 
configuration and between (ls2p3/2) '^-P2 configuration and (Isls) ^5'o configuration with emission of magnetic Ml 
and M2 photons, respectively. By W we denote the transition probability evaluated in this work. The frequency 
of the emitted photon is set equal to ui = Ej — Ep, where Ej, Ep are the energies of the initial and final states, 
respectively. In the case of the nondegenerate levels they comprise the Dirac energies together with one-photon 
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corrections. Accordingly, we do not include the radiative corrections. The contributions of the negative energy states 
to the amplitude are included according to Eqs. (|130p . (jl3ip when performing summation (n) over the entire Dirac 
spectrum. Investigation of the contribution due to the negative energy part of the continuum was performed in 
[27I [5^ . [7^ . For Z > 18 the set g contains configurations in the j-j coupling scheme built on Is, 2s-electrons for 
(1525)^51 configuration and on Is, 2p3/2-electrons for (ls2p3/2) ^-P2 configuration, respectively. For Z < 10, due 
to the poore convergence of the perturbation theory, the set g contains 4000 configurations. Whenever available 
we compare our results with data obtained in other works. The work [isj presents the first relativistic calculation 
of transition probabilities for the (ls2s)'^S'i — > (Isls) transition. The paper by Johnson et al. [50| provides 
a comprehensive review, where the transition probabilities are tabulated for all Z values. However, this work is 
performed neglectin g Q ED effects such as retardation and the contribution arising from the derivative in the vertex 
operator. The work [52l | is performed within the framework of the two-time Green's function method, which is a full 
QED approach, too. However, in this work only the nondegenerate two-electron configurations are considered. The 
present status of measurements of the transition probabilities can be found in [73| . In Table V we give our numerical 
results for transition probabilities between (ls2s) ^Si configuration and the ground state for nuclear charge numbers 
10 < Z < 100. Digits in square brackets indicate the accuracy of the measured values. 

In Tables IHIIIVI we present numerical results for i?l-transition probabilities between (ls2p)'^Pi, ^Pi and ground 
(Isls) ^5*0 two-electron configurations. This provides the first exact QED calculation of the transition probabilities for 
the quasidegenerate configurations. The calculation is performed within the "transverse" and "nontransverse" gauges 
for the emitted photons: Wt and Wm, respectively. The photon propagator was specified in the Coulomb gauge. In 
the case of quasidegenerate levels the energies Ej, Ep include also interelectron interaction corrections of the second 
order. We considered one- and two-photon exchange corrections to the energy levels, taken from (47j . and one photon 
exchange corrections to the eigenvector $„ in Eq. (|146p . The contributions of the negative energy states to the 
amplitude are included in Eq. (I146P in the summation over k (over the complete Dirac spectrum but the set g) and 
in the matrix V (see Eq. (I134p ) in the two-photon exchange corrections. In the columns W^^^ we give the transition 
probabilities calculated with only one photon exchange correction taken into account, i.e., Wn^t is recalculated where 
set Ay(^^ = in Eq. (|135p . The columns W^^P^^ display results of our recalculation of the Wn,t values within 

(2+) 

RMBPT. In the columns AW^ j,' an estimate for the interelectron interaction corrections of higher orders (see the 

end of Subsection IIV Cp is given. The blank fields in the columns Awj^lp' express that the corresponding values are 
smaller than the level of accuracy of the calculation. In the last two columns we giv e the results of Drake [4£] (the 
application of the unified method) and the RMBPT calculations by Johnson et al. [SOl • While in [i^ the "transverse" 
gauge was used, the calculations performed in [soj for (ls2p) ^Pi, ^Pi configurations utilized a "nontransverse" gauge. 
Digits in square brackets again denote powers of 10. 

The diagonalization of the matrix Vu (Eq. (|137p ) implies that we partly take into account the photon exchange 
corrections to all orders. This violates the gauge invariance and, accordingly, explains the deviation between Wn 
and Wt- This deviation also helps us to estimate the contribution of the higher order terms in the expansions Eqs. 
p46p . (|135p . This contribution is larger for small Z values, where the convergence of the perturbation theory in the 
interelectron interaction is poorer. 

The difference between the data in the columns Wn.t and W^^^^'^ determines the nonradiative QED corrections. 
For small Z values the considered configurations are strongly mixed. The transition probabilities for "^Pi levels are 
very sensitive to the mixing matrix B what explains the large value of the QED corrections for the small Z values. 
Note, that the transition probability for the ^Pi level is by several orders of magnitude smaller, than for the ^Pi 
level. This means that the relative correction to the decay of the ^Pi level due to the change of the matrix B is 
essentially larger than the correction to the decay of the level ^Pi. For a large numbers of Z values the mixing of the 
configurations is small and QED corrections appear mainly as QED corrections to (AV2i)fci in the function <!>„ (see 

Eq. isim . 

The perturbation expansion employed for the construction of the matrix V (Eq. (jl34p ) and the one applied in the 
diagonalization of V (Eq. p46p ) are different. The comparison of the columns Wn,t and W^^r^ reveals the importance 
of two-photon corrections to the matrix V for the convergence of the series of the perturbation theory. If we enlarge 
the set g by excited configurations, the values for W^^r^ would approach the ones for Wat.t- However, in order to 
achieve good agreement for small Z, we would have to include about 4000 significant configurations in the set g. The 
matrix V is well investigated for the purpose of evaluating the energies of the configurations [3, [13] , so the evaluation 
of the mixing matrix B in higher orders is a more efficient technique, rather than any enlargement of the set g. The 
technique presented for the calculation of the transition probabilities is a rigorous QED procedure, which allows for 
systematic improvements of the accuracy of the calculation by taking into account corrections of higher orders. 

The accuracy of the presented calculations is determined by the accuracy of the numerical methods, by the contri- 
bution of the omitted orders of the perturbation theories and by the radiative corrections which are not considered 
here. The relative accuracy of the numerical calculation is set to 0.03%. Contribution of the omitted orders of the 
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perturbation theories can be estimated as difference between the values calculated within the different gauges, i.e., 
difference between the and Wt columns. However, this is a very rough estimation, because one of the gauges 
may present better convergence than the other. Partly, the contribution of the omitted orders of the perturbation 

(24-) 

theories can be estimated by the values in columns A j,' ; they show that the "nontransverse" gauge gives consid- 
erably better convergence. To estimate the order of magnitude of the radiative corrections we suppose that they have 
the same order as the other QED effects, i.e., as the difference between the values in Wn,t and W^^^^'^ columns, 
respectively. Accordingly, the values of the transition probabilities calculated within the "nontransverse" gauge (Wn) 
present the most accurate data for the transition probabilities. The estimate of inaccuracy of the data is indicated 
by digits in round brackets. 

Concluding, we can state that at present this paper provides the most extensive and the most accurate calculations of 
the transition probabilities in HCI with intermediate nuclear charge numbers Z. The inclusion of radiative corrections 
into the LPA (which is underway) would yield the most rigorous and powerful approach to the calculation of the 
transition probabilities for HCI with an utmost precision. 
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APPENDIX A: ADIABATIC S-MATRIX 



Adiabatic S-matrix is a modified common S- matrix where adiabatic exponent e"^"'*' is inserted in every vertex. 
The adiabatic parameter Aa is an infinitesimal quantity (Aa — > +0). The presence of the adiabatic exponent switches 
off interaction with the electromagnetic field at t — ±oo. In this appendix we show that the singularities present in 
the common S-matrix (see Section |lTl vanish completely in the adiabatic S-matrix. The singularities arise when one 
makes insertions into the outer electron lines of the Feynman graphs. 

In the present paper we will consider the one-electron ions and the insertions of the self-energy operator. In the 
lowest order of QED perturbation theory S-matrix element corresponding to the process of elastic photon scattering 
on the one-electron ions is given by the Feynman graph in Fig. [TJ We consider the case when electron in the ground 
state oo absorbs photon {k, A), then emits photon {k' , A') and decays back to the ground state. According to energy 
conservation law lo = w' . This graph gives no singularities, accordingly, the adiabatic S-matrix element coinsides with 
the common S-matrix element 



5(0,0) ^ ^(0,0) ^ (_2vrz)5(a>' - co)e' 



(Al) 



The superscripts at the S-matrices indicate the number of insertions of the self-energy operator into the upper and 
lower external electron lines, respectively. Here there are no insertions. 

In the next orders of perturbation theory we have to make insertions of the self-energy operators into the electron 
lines. The insertions into the internal lines yield no singularities. They result in the energy shift of the excited atomic 
states and were investigated in ^45.] . For simplicity of the derivation we omit them. Accordingly, we consider the 
insertions into the outer electron lines. 

The TV insertions of electron self-energy operators into the lower electron line are depicted in Fig. 15. After 
integration over time variables with employment of equality 



-t-oo 



dt 



-Xa\t\+iat 



1 



1 



a + iXa —a + iXa 



(A2) 



the corresponding adiabatic S-matrix element is given by 



S 



iO,N) 



d^Tud^Td^ . . . d^rd2„duJridLL>d^ . . . dudt^dhJsi ■ ■ ■ dus^dfli . . . d^N i'aairu) 



:(-ze)7^"4':''''^*(^»)f E 

ZTT ^ — ' 



'4'niru)tpn{rd^) 
w„ - £„(! - iO) 



x(-ze)7'-^4t')(r,J(z)^ 



27 



1 1 

1 1 



LUn- UJ - UJdi + iXa -^n + UJ + UJd^ + iXa _ 



■ 27r 4^ Wrfi - - »0) 

di 



27r ^ - esi(l - iO) 



1 1 

+ 



1 1 



27r 



1 



1 



ti^dfc - t^Sfc - f^fc + «Aa 
1 



-tJd^ + UJs^ + rifc + «Aa 



1 



CJsfc + f^fc - ti^dfc + i + iAa -Wsfc - ^^fc + Wdfc , 1 + iAa 



d2JV ^^1N-\-\. 



1 



1 



- ^sjv - f^W + ?'Aa. -t^djv + + f^Af + iAa 

1 1 



— £ao + « Aa 



-LiJsN ^ + £ao + *Aa 



(A3) 



In order to make the derivations shorter we will neglect the negative energy part of the Dirac spectrum {i.e., we 
suppose that e„ > 0, £d > 0), since negative energy terms do not generate singularities. 

Consider separately the integral over w-variables and designate it as F. The integrand of F includes all the terms 
of Eq. (jA3P depending on the w- variables: the fractions in the square brackets and the denominators originating from 
the electron propagators. Integration over w„ and ujsi « yields 

2tt 1 



F = dcjdi ■ ■ ■ dujdr. 



1 



1 



Eaa + UJ' - UJ - UJdi + "^iXa Saa + l^' + UJ + UJd^ + 2i\a 

2tt 1 

i {ujdi - £di (1 - iQ)){uJd2 - £si - i^i + iAa) 
1 1 



LOdx - LUd2 + "^iXa -UJdi + ^^d2 + 2iAa 

2tt 1 



i {ujdk - £fe(l - *0))(wd 



k+1 ^Sfc 



1 



1 



2ti 



1 



i {uJdN - £d„(l - jO))(e, 



+ iAa) 
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2l\a 



Integrations in Eq. (|A4p can be performed recursively with the use of equahty 

2tt 1 



1 1 



^dN-l - ^dN + -^dN-1 + ^dN + 2iAa 

27r 1 



1 



f^AT + iX'a) 



27r 



3 



2iX' 



-UJdM + Eao + 2iX' 



1 



1 1 



^ Eqo + 2iA'a + 2iAQ -^djv-l + ^ao + 2iA'a + 2iAa 

1 



(eao - £d„ + 2iA'Q)(£ao ~ - + iX' a) 



where hm -Ra„,a'„ = 0. After integration we get the foUowing expression for Eq. (jA4ll 



2N+1 



I 



I 



1 



(A4) 



(A5) 



(A6) 



(Sao - Edjv + 2jAQ) • • • (Sao - ed„_fc + i + 2fciAQ) • • • (eao ~ ^di + 2NiXa) 

1 



LU + eao-en + {2N+l)iXa 

1 



1 



u - uj' + 2{N + l)iXa -uj + uj' + 2{N + l)iXa 



(A7) 



where lim Rx^ — 0. Function F is singular at when Cd^. — ^ao? k — 1^ . . . ^ N . The term in square brackets in 

Aq - — !-0 



(AS) 



Eq. (|A7P can be written as 



w-w' + 2(A^+l)«Aa -uj + uj' + 2{N + l)iXa 



1 



2(Af+l) 
1 



1 



1 



2(N+1) ^ '-^a 2(7V+1) ^ 

27r 



27r\ / — w' 





2(7V+ 1) V « / V2(A^ + 1) 



(A9) 



Let us restrict ourselves to the case when e^j. = for every k — I, . . . , N (the derivations for the cases when some 
of 7^ can be performed by analogy). Accordingly, we write Eq. (|A7|) as 



2tt 

i 



2N+1 



(ea 



f7jv + iO)---(e 



ao ^SN-fc 



Ojv-fc + «0) • • • (e^o - e^i - fii + «0) 
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2iX„ 



N 



1 

m 

1 



+ (2iV+l)a, 

With employment of Eq. (jA10[) we can write Eq. (jA3|) as 



5 



(CAT) 



(/c',A')* .(fe,A) 



1 



Employing the asymptotic (Aq 



-t-O) equality 



E 

Ar=o 



1 



A 



{x + NiXa)N\ \iXa 



N 



1 



) (^ao) 



N 



2iXn 



■ exp 



A 

iXa 



where 1 2: 1 > lAI we can write 



(AlO) 



(All) 



(A12) 



E 



.(k'X)* .(fe,A) 



X exp 



2iXr, 



(A13) 



(A14) 



Although the condition |a;| > |A| is necessary for Eq. (|A12p . we employ Eq. (|A12p for Eq. (|A11[) for any uj. This is 
considered as an analytical continuation of Eq. (jAll[) to the area close to the resonance and, accordingly, the entire 
complex plane {uj). This analytical continuation was discussed in (l7j . 

If we insert the self-energy operator into the upper (iV„ times) and lower {Nd times) outer electron lines, the similar 
derivations yield 



'nj. ( 



2iXa 



E 



,(fc\A')* .(fc,A) 



1 



) (^ao) 



+ {2N + l)iXa) Ndl 



2iXa 



(A15) 



The value of N in the denominator can be set equal to or Nd without changing the final result since it infiuences 
only the terms which disappear in the asymptotics (Aa +0). Finally, we get 



A„ 



'S'Aa - E 

= {—2TTi)S{ui — uj')e^ 

xE 



.(fc'.A')* .(fc^A) 
^aon ^nao 



X exp 



iXa 



(A16) 
(A17) 

(A18) 
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As the regularized self-energy matrix element for the ground state (ao) has no imaginary part, the absolute value of 
the exponent in Eq. (|A18[) reads 



exp 



) (^ao) 



= 1 



Accordingly, the absolute value of the amplitude is given by 

\U\ = 



,(fc',A')* 4(fe,A) 

in A- ~ P., 



(A19) 



(A20) 



The regularized self-energy matrix element in the denominator is a correction to the energy of the ground state: 

) (^ao)- 



In Eq. (jAlOp we considered only the case when Sd^ = Eoq for every k = 1, . . . ,N . The cases when some of Sd^ ^ 
correspond to the insertions of the second and higher orders self-energy corrections of the "loop-after-loop" type. The 
case when all £d^ 7^ £00 gives the correction to the wave function of the ground state electron cq. 

The presented derivations show that the employment of the adiabatic theory allows for the insertions into the outer 
electron lines within the LPA. They also justify the introduction of the vertex functions $ao, ^"00 (^Aqi ^Aq for 
two-electron ions) in Section Hill The energies of the ground state can be considered to be the full energies, i.e., with 
all the corrections included. 

Finally, we would like to note that the goal of employment of the adiabatic approach was the justification of the 
LPA backgrounds. Formally, the ground state can be investigated within the same matrix formulation of the LPA as 
the excited states which allows to employ the general technique developed in Section HVl 
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TABLE I: Ml-transition probabilities (s ^) between (ls2s)^5'i and (Isls) ^So configurations. The digits in square brackets 
denote the power of 10. 



z 


W 




Ref. [48] Ref. |50] Ref. [52] 


Experiment 


6 


4.867(11) 


1] 


4.856[ 1] 4.860[ 1] 


4.857(11)[ 1]" 


10 


1.097(7) 


4] 


1.087[ 4] 1.092[ 4] 


1.105(18)[ 4]'' 










1.0905(48) [ 4]" 


12 


7.324(30) 


4] 


7.243[ 4] 7.293[ 4] 


7.35(26) [ 4]"* 


16 


1.430(5) 


6] 


1.408[ 6] 1.426[ 6] 


1.422(8) [ 6]" 


18 


4.798(18) 


6] 


4.709[ 6] 4.787[ 6] 




26 


2.078(6) 


8] 


2.002[ 8] 2.075[ 8] 




30 


8.987(21) 


8] 


8.981[ 8] 8.993[ 8] 




50 


1.727(2) 


11] 


1.726[11] 1.729[11] 




54 


3.852(5) 


11] 


3.846[11] 3.856[11] 


3.92(12)[11]-'' 


70 


5.980(10) 


12] 


5.968[12] 5.983[12] 




90 


9.468(23) 


13] 


9.439[13] 9.469[13] 




100 


3.193(1) 


14] 


3.181[14] 





" Schmidt et al. [74]. 
'' Wargelin et al. [75]. 
" Trabert et al. [76]. 

Stefanelli et al. [77]. 
" Crespo Lopez-Urrutia et al. [781 . 
^ Marrus et al. [4]. 



TABLE IL M 2-transition probabilities (s ^) between {ls2p^/2)^ P2 and (Isls) ^5*0 configurations. 

Z W Ref. [50] Ref. [52] 

5 5.016(4)[ 3] 5.014[ 3] 

10 2.258(3) [ 6] 2.257[ 6] 

18 3.145(1)[ 8] 3.141[ 8] 

26 6.515(2)[ 9] 6.510[ 9] 

30 2.104(1)[10] 2.104[10] 2.105[10] 

50 1.365(1)[12] 1.365[12] 1.366[12] 

54 2.560(3) [12] 2.560[12] 

70 2.148(1)[13] 2.146[13] 2.148[13] 

90 1.720(1)[14] 1.718[14] 1.721[14] 

100 4.165(3)[14] 4.156[14] 



TABLE III: _El-transition probabilities (s ^) between {ls2p) ^Pi and (1815)^50 configurations. 



Wn 



AW, 



(2+) 



Wt AVi/j,^+' Ref. [49] Ref. [5Q] 



10 3.096[ 9 
18 1.391[12; 
26 3.925[13; 
30 1.160[14; 
40 6.846[14; 
50 2.104[15; 
60 4.838[15; 
70 9.472[15; 
80 1.672[16; 
92 3.007[16; 



2.917[ 9] 
1.370[12] 
3.898[13] 
1.154[14] 
6.826[14] 
2.101[15] 
4.832[15] 
9.463[15] 
1.671[16] 
3.005[16] 



5.211[ 9; 
1.793[12; 
4.482[13; 
1.258[14 
7.047[14 
2.132[15 
4.874[15 
9.523[15 
1.680[16; 
3.020[16 



4.963[ 9] 5.351(140)[ 9] 0.003[ 9] 5.095[ 9] 0.15[ 9] 5.356[ 9] 5.356[ 9] 

1.772[12] 1.799(6)[12] 1.777[12] 0.011[12] 1.800[12] 1.799[12] 

4.419[13] 4.421(61)[13] 4.396[13] 0.011[13] 4.425[13] 4.421[13] 

1.254[14] 1.251(7)[14] 1.246[14] 0.002[14] 1.252[14] 1.251[14] 

7.041[14] 7.013(34)[14] 6.997[14] 0.007[14] 7.017[14] 7.011[14] 

2.133[15] 2.123(9)[15] 2.120[15] 0.001[15] 2.123[15] 2.120[15] 

4.879[15] 4.855(19) [15] 4.850(15] 0.002[15] 4.853[15] 4.845[15] 

9.538[15] 9.497(26) [15] 9.489[15] 0.003[15] 9.480[15] 9.460[15] 

1.683[16] 1.674(6)[16] 1.673[16] 1.672[16] 1.668[16] 

3.027[16] 3.008(12)[16] 3.006[16] 2.994[16] 
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TABLE IV: Bl-transition probabilities (s ^) between (ls2p) and (Isls) ^5'o configurations. 



z 






p^RMBPT 


^RMBPT 


Wn AW^^'+' 


Wt 




Ref. ^ 


Ref. [50] 


10 


8.607[12] 


8.071[12] 


8.538[12] 


8.107[12] 


8.538(14) [12] 0.007[12] 


8.103[12] 


0.27[12] 


8.851[12] 


8.851[12] 


18 


1.069[14] 


1.052[14] 


1.061[14] 


1.048[14] 


1.061(1)[14] 


1.047[14] 


0.007[14] 


1.071 [14] 


1.070(14] 


26 4.611[14] 


4.578[14] 


4.551[14] 


4.529[14] 


4.553(2) [14] 


4.526[14] 


0.01[14] 


4.570[14] 


4.566(14] 


30 


7.857[14] 


7.815[14] 


7.747[14] 


7. 723 [14] 


7.754(7) [14] 


7. 720 [14] 


0.016[14] 


7.773[14] 


7.763(14] 


40 


2.234[15] 


2.227[15] 


2.211[15] 


2.209[15] 


2.215(4) [15] 


2.210[15] 


0.002[15] 


2.216[15] 


2.212[15] 


50 


5.096[15] 


5.085[15] 


5.064[15] 


5.066[15] 


5.074(10) [15] 


5.065[15] 


0.004[15] 


5.071[15] 


5.057(15] 


60 


1.013[16] 


1.012[16] 


1.009[16] 


1.010[16] 


1.011(2)[16] 


1.010[16] 




1.010[16] 


1.006(16] 


70 


1.819[16] 


1.817[16] 


1.813[16] 


1.816[16] 


1.816(3) [16] 


1.814[16] 




1.813[16] 


1.805(16] 


80 


3.010[16] 


3.007[16] 


3.002[16] 


3.009[16] 


3.008(6) [16] 


3.006[16] 


0.001[16] 


3.000[16] 


2.986(16] 


92 


5.046[16] 


5.043[16] 


5.034[16] 


5.049[16] 


5.045(11)[16] 


5.041[16] 


0.001[16] 




5.001(16] 
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TABLE V: Ml-transition probabilities (s ^) between (ls2s) and (Isls) configurations. The digits in square brackets 
denote the power of 10. 



'7 T^/ 


Z 


W 


Z 


w 




z 


W 




Z W 




21 


2.332(6) [7] 


41 


2.22i[Z) 


[10] 


61 


1.386(2^ 


[12] 


81 2.922(7) 


[13] 




22 


3.757(9) [7] 


42 


2.847(2j 


[10] 


62 


1.644(3) 


r-1 oi 
[12] 


82 3.357(8) 


[13] 




23 5.923(13) [7] 


43 


3.Dz9(2J 


[10] 


63 


1.949(4) 


[12] 


83 3.825(10) 


[13] 




24 9.158(11)[7] 


44 


A a c\c\ ( o\ 


[10] 


64 


2.301(4) 


Fl Ol 

[12] 


84 4.388(12) 


r-i ol 

[13] 




25 


1.391(3) [8] 




0.oU4(oJ 


[lOJ 


65 


2.716(5) 


[12] 


85 5.007(12) 


r-i ol 

[13] 




26 


2.078(6) [8] 


46 


7.283(4J 


[10] 


66 


3.191(6) 


Fl Ol 

[12J 


86 5.679(14) 


[13] 




27 


3.058(6) [8] 


47 


9.102(5) 


[10] 


67 


3.745(7) 


Fl Ol 

[12J 


87 6.499(16) 


[13] 




28 


4.438(7) [8] 


48 


1.132(1) 


[11] 


68 


4.387(8) 


Fl Ol 

[12] 


88 7.350(19) 


Fl ol 

[13] 




29 


6.357(9) [8] 


A n 

49 


1.402(2) 


[11] 


69 5.136(10) 


ri ol 
[12] 


89 8.400(22) 


[13] 


10 1.097(7) [4] 


30 


8.987(21) [8] 


50 


1.727(2) 


[11] 


70 5.980(10) 


Fl Ol 

[12] 


90 9.468(23) 


[13] 


11 2.9DD(loj[4J 


31 


1.259(2) [9] 


51 


2.125(2) 


[11] 


71 6.980(12) 


Fl Ol 

[12] 


91 1.081(3) 


[14] 


12 7.324(30) [4J 
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1.742(2) [9] 


52 


2.600(3) 


[11] 


72 


8.101(15) 


ri ol 
[12] 


92 1.216(3) 


[14] 


13 l.D/o(/)[oJ 


33 


2.388(2) [9] 


53 


3.172(4) 


[11] 


73 9.419(18) 


Fl Ol 

[12] 


93 1.387(4) 


[14] 


14 3.609(13) [5J 


34 


3.244(3) [9] 


54 


3.852(5) 


[11] 


74 


1.089(2) 


Fl ol 

[13] 


94 1.555(4) 


[14] 


15 7.354(23) [5] 


35 


4.367(4) [9] 


55 


4.671(7) 


[11] 


75 


1.260(2) 


[13] 


95 1.774(5) 


[14] 


16 1.430(5) [6] 


36 


5.831(5) [9] 


56 


5.637(8) 


[11] 


76 


1.455(3) 


[13] 


96 2.004(6) 


[14] 


17 2.667(9) [6] 


37 


7.726(6) [9] 


57 6.789(11) 


[11] 


77 


1.677(4) 


[13] 


97 2.262(7) 


[14] 


18 4.798(18) [6] 


38 1.016(1)[10] 


58 8.144(13) 


[11] 


78 


1.931(4) 


[13] 


98 2.550(8) 


[14] 


19 8.358(15) [6] 
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FIG. 1: Feynman graph, describing the photon scattering on an atomic electron. The wavy lines with the arrows describe the 
absorption and emission of photons with momenta k, k' and polarizations e, e', respectively. The double solid line denotes the 
electron in the field of the nucleus, ao corresponds to the ground electron state. 
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FIG. 2: Feynman graph corresponding to the electron self-energy insertion into the electron propagator in Fig. [T] The wavy 
line denotes the virtual photon. The other notations are the same as in Fig. [T] 
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FIG. 3: The process of double-photon scattering on the artificial "lower than ground" state a*. Notations are the same as in 
Fig.[l] 
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FIG. 4: The Feynman graph representing the process of elastic photon scattering on one-electron ion. The double solid line 
denotes the electron in the field of the nucleus (Furry picture of QED). The boxes with the wavy lines describe the absorption 
and emission of the photons by electron the ground state. The letter I at the internal electron line implies the resonance 
approximation, where only the resonant state / remains in the sum over the intermediate states in the electron propagator. 




FIG. 5: The insertion of the one-loop electron self-energy in the internal electron line in Fig. [J] The wavy line describes the 
virtual photon. The other notations are the same as in Fig. U 



ao 



FIG. 6; The Feynman graph, illustrating Eq. (|31|l . The box with the letter E inside corresponds to the infinite number of 
successive insertions of the type Fig. O 
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FIG. 7; The Feynman graph representing the process of the photon emission. The labels / and F correspond to the initial and 
final states. 
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^\A/\J|(d) 

FIG. 8: The Feynman graphs representing the process of photon emission in the LPA. This graph incorporates the graph in 
Fig. [7] The upper, central and down vertices are specified by corresponding subscripts {u), (c) and (d), respectively. 
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FIG. 9: The Feynman graphs representing the single photon transition in a two-electron ion to lowest order in a. The notations 
are the same as in Figs. |4]|8l li {i = 1,2) and Fi (i = f , 2) denote the initial and final states for the two electrons. 



fflO 

k',e' 



fco.eo 



k , e 



Oo 



bo 



bo 



(a) 



^2 



k', e! 



fco, eo 



k , e 



o-o bo 
(b) 



FIG. fO: The Feynman graphs representing the process of elastic photon scattering on the two-electron ion within the LPA. 
These graphs incorporate the graphs in Fig. [9] See in the text the notations ao,6o. 
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FIG. 11: The Feynman graphs representing the process of elastic photon scattering on the two-electron ion. The boxes with 
the letter Aq inside and wavy lines depicts complicated vertices describing the absorption and emission of a photon by the 
two-electron ion. The photon line in the center denotes the emission of a photon with frequency uiq = |feo| corresponding to 
the transition energy from the initial to the final two-electron state. 
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FIG. 12: The Feynman graph representing the process of elastic photon scattering on the two-electron ion. In addition to Fig. 
[TT]the box with letter H inside and with an external photon indicates the complicated vertex for the emission of a photon uio 
corresponding to the transition from the initial to the final state {I ~* F). 
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FIG. 13: The Feynman graph representing the process of elastic photon scattering on the two-electron ion. This graph includes 
the one photon exchange correction and, accordingly, it represents the next order of the perturbation theory (photon exchange 
correction) compared to the graph Fig. 1111 
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FIG. 14: The Feynman graphs representing the process of elastic photon scattering on the two-electron ion. In the lowest order 
of the perturbation theory these graphs reduce to the graphs Fig. 1131 
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FIG. 15: The Feynman graph representing the process of elastic photon scattering on the one-electron ion. Multiple insertions 
of the self-energy operator into the lower outer electron line are made in the framework of the adiabatic approximation. The 
break in the electron lines denotes the possible multiple insertions. 
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